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I.  INTRODUCTION 

1.1  Project  Summary 

The  project  period  covers  IS  August  1992  throu^  31  December  1994.  The  speciHc 
objectives  are  to  investigate  (1)  die  dynamics  of  die  vortex  and  flame  interaction  in  jet  diffusion 
flames  and  (2)  the  transport  phenomena  of  injection  of  supercritical  fluids  into  a  sub-  or  super- 
cfidcal  oivironnMnt  The  overall  objective  of  the  proposed  project  is  to  contribute  to  the  Air  Fcmce 
in  its  effort  of  establishing  a  technology  base  for  fuel  and  combusdcxi  systems  to  meet  the 
requirements  of  future  propulsion  systems.  The  present  research  has  emphasized  on  die  dynamics 
of  the  transport  processes,  for  example,  regarding  the  vortex  and  flame  (or  flow)  interaction  and 
die  injection  of  supercritical  fluids  into  an  environment  which  is  composed  of  dissimilar  fluids  and 
maintained  at  subcritical  states  oi  the  injected  fluid. 

This  report  is  prepared  to  summarize  the  acconqilishmoits  during  the  project  period, 
including  (1)  quantification  of  the  vortex-flame  inter^on  in  a  near-laminar  jet  diffusion  flame,  (2) 
verification  of  the  numerical  simulation  Oi  vortex  structures  inside  a  transititxial  jet  diffusion  flame, 
(3)  experimental  observation  of  supercritical  sprays. 

1.2  Personnel 

The  Project  Director  was  L.  D.  Choi,  Professor  and  Chairperson  of  the  Department  of 

Mechanical  Engineering  at  The  University  ctf  Iowa.  During  die  eariy  stage  of  the  Project,  K.-Y. 

Hsu  served  as  a  post-doctoral  associate  and  conducted  the  experiments  quantifying  the  vortex  and 

flame  interaction  at  Wright  Laboratory.  Hsu  is  pr^ently  a  research  scientist  with  Systems 

Research  Laboratories--  a  Division  of  Arvin/Calspan  and  is  working  at  the  Wright  Laborattny. 

Two  doctoral  students,  S.  Lee  (presently  with  (Cummins  Engine,  Inc.)  and  H.-W.  Lin  (presendy 

with  Matusi-USA)  were  also  funded  during  the  initial  period  of  the  Project,  and  their  doctoral 

degrees  have  been  conferred.  Two  doctoral  students  are  presendy  funded  by  the  Project,  P.  Sui 

who  is  undertaking  measurements  to  quantify  the  supercritical  sprays  and  K.-Z.  Jia  who  is 

w(»king  on  the  trapped  vortex  combustor-  an  innovative  amcept  prc^iosed  by  W.  M.  Roquemore 

of  Wright  Laboratory.  Qose  collaboration  with  the  Wright  Laboratory  has  been  accomplished 
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through  the  technical  management  of  W.  M.  Roquemoie. 

1.3  Liirt  of  Publications  and  Presentations 

During  the  project  poiod,  an  entry  (for  visualization  of  the  vcatex-flame  interaction)  to  the 
1992  Galleiy  of  Fluid  Flow  of  the  Division  of  Ruid  Dynamics,  American  Physical  Society  has 
been  selected  and  published  in  Physics  of  Ruids  A.  (Hsu  et  al.,  1993),  a  manuscript  reporting  the 
inner  vmtex  crossing  frequency  will  be  submitted  to  a  journal  for  publication  (Lee  et  al.,  1994). 
Two  piqrers  were  presented  in  professional  meetings  to  report  the  voitex-flame  interactions  (Hsu  et 
al.,  1993)  and  the  chaotic  aspects  of  the  flame  flicker  in  diffusion  flames  (Sui  et  al.,  1993).  Two 
RiD.  dissertations  were  also  completed.  Lee  (1992)  reported  the  verificatitm  of  the  crossing 
frequoKy  of  die  iiuier  vortices  in  transitional  jet  diffusion  flames  and  Lin  (1993)  reported  the 
effects  of  the  boundary  conditions  to  the  flicker  phenomentxi  of  a  buoyant  jet  diffusion  flame.  A 
list  of  these  publications  foUows. 


ffflimalAaisIffs 

1.  Hsu,  K.-Y.,  Chen,  L.-D.,  Katta,  V.  R.,  Goss,  L.  P.,  and  Roquemore,  W.  M.,  1993, 
"Flame-Vonex  Interacticms  in  a  Driven  IMffesion  Rame,"  APS  Physics  of  Ruids  A.  Ruid 
Dynamics  £in  Gallery  of  Ruid  Motk»),  Vcd.  30,  p.  S4. 

2.  Lee,  S.,  Chen,  L.-D.  and  Roquemore,  W.  M.,  1^4,  "Vortex  Crossing  Frequency  in  a 
Transitional  Jet  Diffusion  Rame,"  to  be  submitted  to  Combustion  and  Rame  (1994). 

Conference/Proceedings 

1 .  Hsu,  K.  Y.,  Chen,  L.-D.,  Katta,  V.R.,  Goss,  L.P.,  and  Roquemore,  W.  M.,  1993, 
"Ej^oimental  and  Numerical  Investigations  of  die  Vcvtex-Flame  InteracticMis  in  a  IMven 
Jet  Diffusion  Flame,"  AIAA  Pqier  93-0455,  Thirty-first  Aerospace  Science  Meeting  heW  in 
Reno,  Nevada,  11-14  January,  1993. 

2.  Chen,  L.-D.,  Sui,  P.C.,  and  Roquemore,  W.M.,  1993,  "Dynamic  (Characteristics  of  a 
Buoyant  Jet  Diffusion  Flame,"  Bulletin  of  the  American  Physical  Society.  Vol.  38,  No.  10. 

Thesis 

1 .  Lee,  S.,  "Row-Flame  Interaction  of  Nitrogen-Diluted  Hydrogen  Jet  Diffusion  Rame," 
RlD.  Dissertation,  Department  of  Mechanical  Engineering,  The  University  of  Iowa,  Iowa 
City,  Iowa,  December  1992. 

2.  Lin,  H.-W.,  "Numerical  Simulation  of  the  Dynamics  and  InstaNlity  of  Rame  Flicker  when 
Subject  to  Perturbed  Boundary  Conditions,"  Ph.D.  Dissertation,  Department  of  Mechanical 
Engineering,  The  University  of  Iowa,  Iowa  City,  Iowa,  December  1993. 
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II.  VORTEX-FLAME  INTERACTION 

The  initial  results  of  the  vortex-flame  interaction  study  are  reported  in  a  recent  paper 
presented  at  the  thirty-first  Aerospace  Science  Meeting  held  in  Reno,  Nevada,  1 1-14  January,  1993 
(Hsu  et  al.,  1993).  A  specific  type  of  the  vortex-flame  interaction  is  studied,  in  that  the 
characteristic  length  of  the  vortex  is  much  larger  than  the  luminous  flame  thickness  and  the  vortex 
characteristic  time  is  much  longer  than  the  reaction  characteristic  time.  This  vortex  was  femned  due 
to  an  acoustic  driving  to  the  fuel  jet  (in  a  co-flowing  air  envirraiment)  at  30  Hz.  The  acoustic 
driving  periodically  reproduces  a  large  ring  vortex,  and  an  accompanying  counter-rotating  ring- 
vtntex  pair.  These  large  ring  vortices  interact  with  tte  thin  flame  typical  of  a  near-laminar  jet 
diffusion  flame.  This  interaction  leads  to  local  extinction  through  stretched  flame  characteristics. 
The  preliminaiy  data  obtained  thus  far  suggest  that  partial  premixing  may  have  occurred  inside  die 
ring  vortex  pair.  It  also  suggests  that  convoluted  flame  sheets  or  distributed  reactions  may  exist  in 
the  vortex  pair.  The  above  speculation  needs  to  be  confirmed  by  e}q)eriments.  The  preprint  of  the 
piqper  by  Hsu,  et  al.  (1993)  is  included  in  Appendix  A. 
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III.  NUMERICAL  INVESTIGATION  OF  A  JET  DIFFUSION  FLAME 

A  manuscript  (Lee  et  al.,  1994)  quantifying  the  crossing  frequency  of  the  inner  vntices  and 
verification  of  the  numerical  simulation  has  been  prepared  and  will  be  submitted  for  publication  in 
an  archival  journal  The  crossing  frequency  of  the  inner  vortex  in  a  transitional  jet  diffusion  flame 
was  measured  as  a  function  of  axial  location  using  a  line-visualization  technique.  Step-wise 
decreases  oi  the  crossing  frequency  were  observed  in  the  experiments,  in  qualitative  agreement 
with  that  reported  in  the  literature.  The  numerical  simulation  based  on  a  simple  flame-sheet  iiKxiei 
was  shown  to  reproduce  the  inner  vortex  in  the  transititmal  jet  diffusion  flame  studied.  Random 
perturbations  applied  at  the  fiiel-jet  exit  were  required  for  the  prediction  of  roll-up  vortices  inside 
die  flame.  The  simulation  results  also  revealed  a  step-wise  decrease  of  the  crossing  frequency 
along  the  axial  ditectitxi.  The  simulation,  although  ctqitured  die  major  ctxnponents  of  the  crossing 
frequency  in  the  jet  near  field,  was  only  in  qualitative  agreement  widi  the  experimental  crossing 
frequency  data  as  a  function  of  axial  location.  The  Strouhal  number  scaling  of  the  vortex  crossing 
frequency  showed  that  the  Strouhal  number  was  bounded  by  0.2  and  1.0,  typical  of  the  Kelvin- 
Helmholtz  instability  of  cold  jets.  The  results  obtained  during  the  project  period  suggested  that  the 
inner  vortex  is  a  Kelvin-Helmholtz  type  of  the  convective  instabili^  which  is  modified  due  to  the 
inesoKe  of  the  high-tenqierature  flame  zone.  The  two-dimensional  numerical  model,  di^fcae,  is 
capable  of  rqmxiucing  the  frequency  characteristics  of  the  buoyancy-induced  vortex  in  transitional 
jet  diffusion  flames  as  well  as  capturing  die  frequency  con^xments  of  a  Kelvin-Helmholtz  type  of 
instability  inside  the  flame.  A  cqiy  of  the  manuscript  t^Lee  et  al.  (1994)  is  included  in  Appendix 
B. 


IV.  SUPERCRITICAL  SPRAYS 
4.1  Experimental  Considerations 
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The  injection  high-ixessuie  sprays  (i.e.,  die  injected  fluid  pressure  being  around  or 
higher  tiian  its  themxxlynamic  critical  pressure)  can  be  cataloged,  for  example,  into  three  regimes; 
high  pressure  sprays  (Regime  1),  supercritical  sprays  (Regime  2)  and  supercritical  flashing 
injection  (Regime  3),  shown  in  the  tenqterature-entropy  diagram  (or  T-s  diagram)  of  Hg.  1.  If  one 
cmisiders  the  spray  into  an  environment  being  con^sed  of  its  own  viyx>r,  the  T-s  diagram  shown 
in  Fig.  1  adequately  describes  the  spray  regimes.  Regime  1  signifies  such  a  typical  fuel  injection 
condition  in  a  diesd  engine  or  in  a  conventional  jet  engine.  The  fuel  at  the  injector  exit  usually  has 
a  supercritical  pressure  but  a  sub-critical  temperature.  Therefore,  the  atomization  and  vaporization 
is  crucial  to  the  mixing  of  this  injection  process.  In  Regime  2,  the  injector  fluid  is  typically  super¬ 
critical  (but  near-critical).  For  a  single^on^xment  system,  the  mixing  of  doise  fluids  dictates  the 
tranqxxt  process.  However,  when  a  dissimilar  fluid  or  a  multi-component  system  (such  as  the 
combustor  condition)  is  involved,  the  thermodynamic  state  varies  as  the  local  mixing  omdition 
varies.  For  exanqrle,  the  local  thenrKxiynamic  critical  pressure  and  temperature  are  to  be 
determined  by  the  local  multi-conqxxient  phase  equilibrium  when  the  injector  fluid  is  in  contact 
with  a  dissimilar  fluid.  It  is,  therefore,  conceivatde  tiiat  the  mixing  of  a  supercritical  fluid  with  a 
dis-similar  fluid  can  result  in  recondensation  of  the  injects  fluid  when  the  local  phase  equilibrium  is 
maintained.  The  transpot  of  the  injector  supercritical  fluid,  however,  should  determine  die  local 
mixture  states.  The  injection  process  in  Regime  2  yields  a  mixture  diermodynamic  state  located  in 
the  liquid  side  of  the  two-phase  done.  The  injection  process  of  a  supercritical  fluid,  however,  can 
also  result  in  a  flashing  injection  whoi  the  mixture  thermodynamic  state  is  located  in  the  vapor  side 
of  the  two-phase  dome,  for  exarrqile.  Regime  3  of  Hg.  1.  The  flashing  injection  can  be  achieved 
when  the  injector  fluid  has  an  excess  enthalpy  (o  having  a  higher  temperature)  exceeding  the 
enthalpy  of  vaporization  of  the  mixture.  This  concept  was  explored  by  Szetela  and  TeVeltte  (1982) 
and  termed  external  fuel  vaporization  for  aircraft  gas  turtnnes.  The  injection  of  a  supercritical  fluid 
into  a  dissimilar  fluid  environment  is  the  focus  of  the  present  study. 


Figure  1.  A  Schematic  of  High-Pressure  Spray  Regimes 
A  simulant  fluid,  SF6  (Madieson  SCF  purity  grade),  is  ennployed  to  study  the  injection 
process.  This  fluid  is  non-toxic,  nonflammable,  and  environmentally  friendly,  and  it  also  has 
properties  similar  to  ^ical  hydrocarbon  fuels.  Some  noted  properties  are  its  molecular  weight, 
the  near  ctidcal-state  dynamic  viscosity  and  density.  Thoefore,  diis  simulant  fluid  makes  a  good 
candidate  fluid  to  study  the  injection  process  of  liquid  fuels.  The  critical  temperature  and  pressure 
of  SF6  are,  respectively,  318.7  K  and  37.1  bar,  and  those  of  typical  hydrocarbon  fuels  are  700  K 
and  20  bars.  Therefore,  the  heating  requirement  to  bring  SF6  to  a  supercritical  state  is  less 
demanding  than  that  to  bring  hydrocarbon  fuels  to  a  siinilar  thamodynamic  state. 

The  experimental  set-up  includes  a  fuel  delivery  system,  a  spray  nozzle,  and  a  high- 
pressure  spray  chamber  (e.g.,  see  Fig.  2).  A  HFLC  pump  (Gilson  305  Piston  Pump,  with  pump 
head  Model  50SC)  is  used  to  deliver  the  injector  fluid  at  preset  flow  rates.  The  pump  pressure, 
however,  can  not  be  accurately  pre-set  because  it  is  determined  by  the  pressure  drop  needed  to 
deliver  the  desired  flow  rate  at  the  chamber  pressure.  The  high-pressure  chambo'  (American 
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Research  CcHporation)  has  an  inside  diameter  of  102  mm  and  height  of 430  imn.  The  total  volume 
oi  the  chamber  is  i^proximately  3.6x10*3  m3.  Four  gauge  glasses  are  mounted  on  the  chamber, 
each  separates  by  90  degrees.  The  high  pressure  chamber  is  designed  for  operation  with  a 
pressure  as  high  as  68  bars  (Roquemore,  1992).  Prior  to  the  injector  exit,  the  SF6  is  heated  by  a 
wax  bath,  with  a  heating  uq)e  wraiq)ed  around  the  stainless  steel  tubing.  A  near-critical 
tenq)erature  is  maintained  at  the  injector  exit  as  can  be^n  from  the  reduced  tenq)aature  (Tr) 
shown  in  Table  1.  Two  spray  injectors  have  been  used,  one  is  a  commercially  available  pressure 
atomizer  (Spraying  System  17S30-SS(X)S)  with  a  specified  orifice  size  (127  pm  in  diameter)  and 
the  other  was  a  Wright  Laboratory  designed  and  tiuuiufactured  glass  injector  (Switzer,  1992).  The 
glass  nozzle  has  a  long  passage  of  1/di  equal  to  36  (1  being  the  passage  length  and  di  the  inside 
diameter)  and  a  sudden  contraction  prim  to  die  nozzle  exit;  the  contraction  ratio  di/de  is  equal  to  6, 
de  being  the  injectm  diameter.  The  diameter  of  this  nozzle  is  estimated  to  be  150  pm.  Thenozzle 
is  installed  vertically  in  the  top  center  of  the  chambm.  A  sub-miniature  k-type  thermocouple 
(Omega  KMQSS-020E)  is  inserted  into  the  nozzle  assembly  at  about  20  mm  upstream  of  the  nozzle 
exit 

The  SF6  was  injected  into  the  spray  chamber  fiUed  with  N2  (designated  as  SF6/N2),  CO2 
(designated  as  SF^/CXDj),  m  its  own  vapor  (desi^ated  as  SF6/SF6).  The  SFe/SFe  condition 
serves  as  a  reference  for  the  flashing  injection  of  a  supercritical  fluid.  The  mass  flow  rate  of  SF6  is 
maintained  constant,  3.75  x  10*^  kg/s,  through  the  tests.  It  should  be  noted  that  there  is  no 
throttling  process  in  the  flow  passage  after  the  pump.  Wax  bath  and/or  heating  tape,  or  nichrome 
wire  are  used  to  heat  die  SF6.  The  test  chamber  is  charged  by  CO2.  N2,  or  SF5  vapor  as  the 
ambient  gas  to  desired  pressures  prior  to  the  experiments.  A  small  amount  of  the  ambient  gas  is 
allowed  to  flow  through  the  chamber  during  the  tests.  The  chamber  pressure  is  maintained 
constant  and  is  contrcdled  by  adjusting  the  fine  metering  valve  at  the  bottcxn  of  the  chamber.  The 
test  conditions  are  summarized  in  Table  1.  Hie  test  conditions  are  determined  by  maintaining  a 
same  volumetric  flow  rate  (mass  flow  rate)  fn*  the  injection  conditions.  Therefore,  a  direct 
compariscm  of  drc^let  formation  (reccmdensation)  and  spray  length  can  be  made.  The  volumetric 
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flow  rate  is  determined  to  yield  a  turbulent  Reynolds  number  (Re)  at  the  injeaw  exit  and  a  spray  in 
the  atcxnization  regime  following  the  jet  breakup  regime  boundaries,  for  example,  summarized  and 
defined  by  Reitz  and  Bracco  (1986).  The  spray  regime  is  identified  by  the  Ohnescxge  (Oh)  and 
Reynolds  numbers  based  on  the  injector  exit  condition.  The  Froude  number  (Fr)  of  the  test 
conditicxis  is  in  the  range  7.3  x  10’'^  to  3.3  x  10“^,  suggesting  that  these  conditions  are  momentum 
dominated.  The  chamber  is  maintained  with  a  near-critical  pressure  as  can  be  seen  from  the 
dimensionless  pressure  (Pr'),  which  is  normalized  by  die  critical  pressure  of  SF6  shown  in  Table 
1.  The  dimensionless  pressure  of  the  chamber,  Pr',  is  also  the  reduced  pressure  of  SF6  at  the 
injector  exit  It  should  be  menticxied  that  no  heating  is  applied  to  the  chamber.  Therefore,  the 
chamber  is  assumed  to  be  at  ixxxn  temperature  (or  298  K)  for  all  the  conditions,  or  having  a 
corresponding  reduced  temperature  based  on  the  SF6  critical  temperature  (Tr*)  of  0.92. 

Two  diagnostic  techniques  are  employed,  "planar"  imaging  and  shadowgraph  visualization. 

The  "planar"  imaging  was  accomplished  by  passing  through  the  center  of  the  spray  with  a  vertical 

laser  sheet  A  He-Ne  laser  (10  mW  at  632  nm)  is  used  as  die  illuminating  laser  source.  The  spray 

image  is  recorded  by  a  CCD  video  camera  (Cannon  A-1,  or  Panasonic  AG-180U)  which  is 

positioned  normal  to  the  illuminating  light  sheet  The  video  image  is  digitized  by  a  frame  grabber 

(Scion  Image  Capture  11;  8-bit  resoluticxi)  interfaced  to  a  micro-computer  (Mac  IIx)  for  analysis 

» 

and  measuremoits  of  the  spray  length.  Limited  experimrats  are  conducted  using  a  pulsed  (10  ns) 
Nd:YAG  laser  (wavelength  at  532  nm)  as  the  illuminating  light  source  to  resolve  the  instantaneous 
flow  features.  The  NdrYAG  lasn  is  used  in  conjunction  with  an  intensified  CCD  (ICCD)  camera 
positicMied  normal  to  the  illuminating  light  sheet  The  ICCD  camera  is  based  on  a  commercial  unit 
(PhotcHnetrics  Model  CH220/camera  head  and  Thomson  TH7882  CCD  chip,  384  x  576  pixels) 
with  a  fiber-optically  bonded  intensifier  (ITT  Model  F41 1 1  S-20  intcnsifier). 

The  thickness  of  the  illuminating  laser  sheet  is  about  250  p.m  in  thickness.  TTiis  thickness 
can  not  provide  spatially-resolved  results  of  the  near-fieki  imaging  because  of  being  on  the  similar 
order  of  magnitude  of  die  injector  diameter  (around  130  pm).  Therefore,  "light  scattering"  instead 
of  "planar"  imaging  will  be  used  in  the  discussion.  To  assist  interpretation  of  the  "light  scattering" 


Table  1.  Summary  of  Experimental  Condition 


d(nm) 

127  (Delevan  Nozzle);  150  (WL  Nozzle) 

Vol.  Flowrate  (m^/s) 

2.5  X  10-7 

2.5  X  10-7 

Mass  Flowrate  (kg/s) 

3.75  X  10  -4 

3.75  X  10  ^ 

Tr 

0.92  -  1.01 

0.92  -  1.01 

0.92  -  1.01 

Pr* 

0.9  -  1.1 

0.9  -  1.1 

0.9  -  1.1 

Re 

29000  -  155000 

29000-155000 

We 

23500  -  oo 

23500 -oo 

Oh  (Ohnesorge  No.) 

5.1  X  10-3.00 

5.1  X  10-3 .  oo 

Fr 

2  m2 

FN  FN  = - 

V  P  ApJ 

(Flow  Number) 

7.3  X  10-7  -  3.3  X  10-6 

2.5  X  10-8 . 6.5  X  10-8 

7.3  X  10-7  -  3.3  X  10-6 

results,  shadowgraph  images  are  taken  using  both  the  He-Ne  laser  as  the  light  source. 
4.2  Results  and  Discussion 
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The  light  scattering  image  is  used  to  determine  the  "spray  penetration  length"  (cm:  the  spray 
length).  The  video  image  is  digitized  by  a  frame  grabber  (8-bit  dynamic  range).  Ten  randomly 
selected  frames  are  used  to  obtain  averaged  intensities.  Typical  axial  intensity  profiles  are  shown 
in  Rg.  3  for  ?F6  injection  into  N2  with  an  injector  exit  temperature  of  Tr  =  0.95, 0.96, 0.97, 0.98, 
0.99  or  l.(X).  The  chamber  is  maintained  at  Pr*  =  1.0.  The  normalized  light  intensity  shown  in 
Fig.  3  is  based  on  an  8-bit  dynamic  range,  i.e.,  255  being  the  highest  intensity  and  0  the  lowest 
The  measured  intensity  decreases  as  the  axial  distance  is  increased,  a  result  of  the  decreasing  light 
scattering  cross-section  due  to  the  decrease  of  droplet  size  or  drc^let  number  density  (or  both). 

Based  cmi  a  threshold  value  5,  the  measured  spray  length  plotted  in  1/d  G  being  the  spray 
length  and  d  the  nozzle  diameter)  is  shown  in  Fig.  4.  The  spray  length  decreases  with  increasing 
Tr  fcMT  SF6  injection  into  CO2  and  N2  environments.  Three  chamber  pressures,  Pr'  =  0.9, 1.0  and 
1.1,  are  examined  in  the  experiments.  The  spray  length  is  shorter  for  lower  chamber  pressures, 
particulaily  true  for  the  injection  into  N2  environmott 


Figure  3.  Axial  Profiles  of  Scattered  Light  Intensity. 


To  assist  inteipretation  of  the  measured  sjnay  length,  a  parabolic  flow  calculation  based  on 
the  locally  hcxnogeneous  flow  (LHF)  model  assumption  in  conjuncdcMi  with  a  k-e-g  turbulence 
model  is  performed.  The  SF($  density  as  a  function  of  Tr  is  shown  in  Hg.  S  for  three  pressures,  Pr 
0.9, 1.0,  and  1.1.  The  density  decreases  rapidly  when  the  critical  temperature  is  approached. 
The  calculated  potential  core  length  (Ip)  using  a  GENMIX-based  c(Hiq)uter  program  is  shown  in 
Hg.  6  for  injection  of  near  and  supercritical  SFg  into  CO2  and  N2  environments  at  conditions 
similar  to  the  experiments.  The  predicted  potential  core,  shown  in  Fig.  6,  decreases  with 
increasing  ten:q)erature.  When  temperature  is  increased,  the  density  of  SF6  decreases  so  does  the 
doisity  ratio  of  SF6  to  N2  or  CO2.  It  is  this  decrease  in  density  ratio  that  yields  a  shorter  potential 
core.  The  measured  spray  penetration  length  (Hg.  4)  shows  a  decreasing  trend  similar  to  diat  of 
the  predicted  potential  ctne  (Fig.  6).  This  comparis(»  suggests  that  the  mixing  of  dense  fluids 
should  be  investigated  in  the  study  of  injection  of  suporritical  fluids. 


0.92  0.94  0.96  0.98  1.00  1.02  1.04 

Tr 


Figures.  Denaty  of  SFg  versus  Tr 
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To  investigate  die  effects  density  ratio  (y)  on  spray  penetration,  data  are  plotted  as  a 

function  of  yin  Fig.  7.  The  critical  tra^ierature  conditions  are  also  marked  in  die  figure  for  the 
dueediamber  pressures  examined.  The  spray  length  decreases  with  decreasing  densiQr  ratio.  For 
the  N2  environment,  a  sharp  decrease  is  observed  at  y  »  25, 18,  and  9,  or  the  criticai  temperature 
state  (Tr  =  1.0)  fcH*  Pr  =  1.1, 1.0, 0.9,  respectively.  This  correlation  with  the  oitical  temperature, 
however,  is  not  observed  for  the  CO2  environment 

The  mixture  pseudo-critical  pressures  and  tonperatures  are  plotted  as  a  function  of  mixture 
fraction,  f.  in  Figs.  8  and  9,  respectively.  The  mixture  fraction  is  defined  as  the  mass  fraction  of 
injected  fluid,  SF6,  in  the  flow  field  with  f = 0  denoting  the  ambient  fluid  and  f  =  1  the  injected 
fluid.  Also  plotted  in  Fig.  9  are  the  mixture  tenqieratures.  An  adiabatic  mixing  between  SF6  and 
ambient  fluid  is  assumed  with  ambient  fluid  at  290  K  and  SF6  at  Tr  =  .92,  .94,  .96,  .98, 1.00,  or 
1.02.  The  tenqierature  is  obtained  based  on  the  Lee-Kesler  equation  of  state  and  the  pseudo-critical 
properties.  Ftom  a  single  coriqionent  pdnt  of  view,  die  mixture  will  be  at  a  supercritical  state 
when  T  >  Tr  and  P  >  Pr.  The  pseudo-critical  pressure  of  the  SF6/N2  mixture  is  lower  than  die  SF6 
critical  pressure  (or  chamber  pressure)  over  the  region  f<  0.6  and  higher  over0.6^f<  1.  For 
SF6/N2,  the  mixture  tnrqierature  is  generally  above  the  mixture  critical  tenqierature,  except  for  f  > 
0.92  when  die  injector  fluid  is  at  Tr  <  1.00.  Thus,  SF6/N2  will  be  generally  in  the  supercritical 
regime  fw  f  <  0.6  and  in  the  subcritical  regime  fix'  0.6  ^  f  <  1.0  (except  for  Tr  <  1  and  f  >  0.92). 
The  pseudo-critical  pressure  of  the  SF6/CO2  mixture  is  higher  than  the  chamber  pressure  for  the 
entue  mixture  fiaction  space  (i.e.,  all  the  mixing  states).  Therefore,  SF6/CO2  will  be  in  the 
subcritical  regime  through  out  the  mixing.  It  is  also  wmth  noting  that  the  SF6^G02  mixture 
traqierature  is  ^nerally  lower  than  the  pseudo-critical  temperature,  except  for  mixture  fiactions 
near  unity  when  the  injector  fluid  is  at  Tr  >  1  (i.c.,  for  f  >  0.94  with  injector  fluid  at  Tr  =1.02 
among  the  six  injector  temperatures  examined). 

The  mixture  pseudo-critical  prqierty  calculation  suggests  that  die  subcritical  state  exists  in 
SF5/N2  quays  over  the  mixture  fiaction  bounded  by  0.6  ^  f  <  1  and  exists  in  SF6/CX)2  sprays 
over  the  entirely  mixture  fraction  space  except  for  die  small  regime  near  f  =  1  as  noted  earlier.  The 


Teiiq>erature  (K) 


Hguie  9.  Calculated  Mixture  Pseudo-Critical  Teiuperatures  as  a  Function  of  Mixture  Ruction. 
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SF6/C02 


Tr*1.01  Tr*1.00  Ti^.99  Tn=0.93 
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Figure  11.  Instantaneous  images  of  SF6/CC)2  sprays  at  Pr  =  1.0 
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SF6  Injection  into  N2  @  Pr=  1.0 

Shadowgraph  Sheet>lighting  Shadowgraph 

Tr=0.94  Tr=0.94  Tr=1.00 


Sheet-lighting 

Tr=1.00 


Figure  12.  Instantaneous  Image  of  SF6/N2  Sprays  at  Pr  =  1.0. 
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calculation  tlwrefwe  suggests  that  the  SF5^CX)2  spray  has  a  wider  spray  angle  and  a  longer  spray 
length  under  a  comparable  injector  condition.  The  wider  spray  angle  and  longer  spray  length 
indeed  oui  be  seen  from  the  avoaged  light  scattering  image  shown  in  Hg.  10.  When  the  spray 
length  is  plotted  as  a  function  of  density  rado  as  shown  by  Fig.  7,  the  CO2  environment  with  Pr'  = 
1.0  (i.e.,  with  die  chambo'  being  maintained  at  SFg  cridcal  pressure)  has  a  longer  sjnray  length 
ovw  12  <  Y  <  18.  The  CO2  environment,  however,  has  a  shorter  spray  length  fw  Pr*  =  0.9  over 
10  <  Y<  18,  in  contrary  to  die  predicdon  of  pseudo-cridcal  mixture  states.  Only  when  the  injected 
fluid  density  decreases  below  the  cridcal-sate  value  (i.e.,  y  below  9  fw  Pr*  =  0.9)  the  N2 
environment  has  a  shorter  spray  length.  At  diis  density  rado  (or  rather  when  the  injector  fluid 
tenqierature  reaches  the  cridcal  tenqierature),  the  spray  length  decreases  sharply.  A  similar  sharp- 
decrease  is  also  seen  fcv  the  N2  environment  at  Pr'  =  1,  e.g.,  at  y**  18  in  Fig.  7.  As  for  CO2 
envircMiment,  such  sharp  decrease  is  less  obvious.  It  repeats  that  die  injecdon  of  near-  or  siqier- 
cridcal  SF6  into  N2  environment  results  in  a  flashing-lilce  atomizadon  while  the  CX^  environment 
does  not 

The  instantaneous  light  scattering  images  of  the  SF5  qvays  in  CC^  environment  are  shown 
in  Hgure  11  ftv  injector  fluid  at  Tr  =  0.93, 0.99, 1.00  or  1.02.  The  calculadon  of  mixture  pseudo- 
critical  stales  of  SF6fC02  sprays  suggests  diat  the  mixture  fiacdon  is  typically  subciidcaL  The 
spay  imaging  of  Fig.  11  shows  jet-like  mixing  or  strwtures  ftv  die  CX32/SF6  sprays.  The  jet-like 
structures  are  shown  by  die  "fine  mists"  formed  in  the  flow,  in  which  large  vordcal  structures  can 
also  be  identified.  When  the  injector  tonperature  increases,  die  "mist"  region  is  reduced  which  is 
evidoit  from  die  comparison  ci  the  images  taken  at  Tr  =  0.93  and  1.02.  The  formadon  of  mist-like 
droplets  in  die  flow  could  result  from  the  local  equilitmum  and  it  in  turn  yields  suboidcal  state- 
dn^.  It  appears  that  the  mixing  process,  instead  of  the  break-up  of  a  liquid  core,  dictates  the 
fcnmatkxi  and  die  "q;iread"  of  the  drcqilets.  When  the  injector  temperature  is  increased,  the  injected 
fluid  aldiough  inidalfy  being  siqieicridcal,  the  mixture  is  sub-<ridcal.  Therefore,  droplets  will  f(»m 
when  die  mixing  between  die  injector  and  ambient  fluids  occurs.  Far  away  from  the  injectmr,  e.g., 
f  s  0,  the  ambient  fluid  exists  in  viqior  phase.  Moving  toward  the  jet  core  (or  f  >  0),  the  liquid 
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{diase  (or  fonnation  of  dn^lets)  will  i^pear  when  the  local  SF^  vapw  pressure  (or  fugacity) 
assumes  the  value  determined  by  the  phase  equilibrium  and  mixing  considerations.  When  the 
injector  teiiq)erature  (or  the  excess  enthalpy)  increases,  the  spray  bdiaves  like  a  flashing 
atomization  df  supo^eated  liquids.  When  the  temperature  or  the  superheat  increases,  the  available 
excess  enthalpy  for  viqx)rizati«i  increases  and  the  volume  occupied  by  the  liquid  phase  decreases. 

As  for  the  injection  of  SF6  into  N2  environment,  the  mixture  fraction  is  supercritical  for  f  < 
0.6  and  strcmg  flashing-like  atomization  is  likely  to  exist  when  the  injector  fluid  is  maintained  at  a 
near  or  supoxnitical  tenq)erature  (cf.  Figs.  8  and  9).  The  instantaneous  imaging  of  SF6/N2  sprays 
shows  (Fig.  12)  diat  the  liquid  phase  (or  a  stronger  scatter)  exists  only  in  the  core  of  the  jet  The 
mist-like  droplets  are  not  obvious  in  die  ICCD  image  of  Fig.  12;  however,  the  jet  boundary 
extends  Off  beyotKl  the  dense  erne  as  evident  from  the  shadowgraph  images.  The  absence  of  mist- 
like  droplets  is  probaUy  due  to  the  fact  that  the  mists  or  drc^lets  are  indeed  not  formed  or  the 
intrasity  of  scattered  li^t  is  below  the  detecticHi  limit  set  for  the  particular  visualization.  When  the 
injector  ten^erature  is  increased  from  Tr  =  0.94  to  l.(X),  the  volume  of  the  dense  core  decreases 
(shown  by  Ught  scattering  images)  and  the  vapor  core  remains  qualitatively  similar. 

4.3  Uncertainty  Estimates 

The  sensitivity  in  die  measurements  of  the  spray  length  due  to  the  threshold  value  is 
summarized  in  Hg.  13.  The  relative  error  is  defined  as  the  length  for  pixel  intensity  below  a 
dueshedd  value  divided  by  the  measured  length.  The  threshold  lengdi  is  determined  by  a  cubic 
polynomial  curve  fit  of  die  axial  profile  and  the  threshold  pixel  value  is  chosen  to  be  5  (the  full 
scale  value  is  255).  Most  errors  are  within  5%  of  the  measured  lengdi.  The  errcff  becomes  largo* 
udioi  die  spray  length  is  short  due  to  strong  flashing-like  atomization,  e.g.,  SF6/CO2  at  Pr  =  0.9, 
or  when  the  dense  fluid  is  confined  in  the  core  regicxi  and  the  break-up  of  fluid  parcels  is  recorded, 
e.g.  SF6/N2. 
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Tr 

Figure  13.  Entn*  Estimate  for  the  Spray  Length  Measurements. 


4.4  Summary 

The  spray  length  measurement,  the  calculation  of  mixture  pseudo-critical  temperatures  and 
pressures,  and  the  instantaneous  light  scattering  and  shadowgraph  images  provide  new  insight  into 
the  injecticMi  of  supocritical  fluids  into  an  environment  of  dissimilar  fluids.  A  simulant  fluid,  SF6, 
and  two  envinximents,  N2  and  CX)2>  are  examined.  Mist-like  droplets  with  characteristics  of  jet 
mixing  are  ^ical  for  the  CX)2  envinmment,  ccmsistent  with  the  mixture  critical  state  calculatirai 
which  shows  that  the  respective  mixture  fraction  space  is  subcritical.  Narrow  and  dense  cores  are 
typical  in  the  light  scattering  imaging  of  SF^  in  N2  environment  with  the  shadowgraph  showing  a 
jet-like  spread,  consistent  with  the  mixture  critical  state  calculation  showing  that  the  mixture 
firiKrtion  space  is  subcritical  primarily  fra-  the  region  0.6  ^  f  <  1.  The  mixing  of  dense  fluids 
tqrpears  to  dictate  the  spread  of  the  mist-like  drcq)lets  in  CO2  environment  and  the  spread  of  the 
vapor  in  N2  environment  The  flashing-like  atomization  is  another  characteristics  in  the  sprays 
reported  in  the  present  investigaticm.  The  present  results  call  for  further  investigation  on  the 
mixing  characteristics  of  supercritical  sprays. 
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Abstract 

This  paper  presents  the  initial  results  of  a 
vortex-flame  interaction  that  is  periodically 
refU'oduced.  The  interaction  has  features  that  are 
similar  to  those  observed  in  a  turbukat  jet  diffusion 
flame.  An  acoustic  speaker  is  us^  to  form 
axisymmetric.  counter  rotating  ring  vortices  inside 
the  flame  surface.  Phase-locked  Reactive-Mie- 
Scattering  (RMS)  and  Planar-Laaer-Induced- 
Fluoiescence  (PLIF)  ate  used  to  visualiae  the  vortex 
structure  and  OH  zone,  respectively.  The  evoludon 
of  the  vortex-flame  interaction  is  presemed.  The 
forcing  vrmex,  convected  in  the  radial  direction, 
causes  the  flame  to  bulge  outward.  When  the  flame 
surface  is  further  stretched,  a  local  quenching  around 
the  flame  is  observed.  Thin-Filaroeat-Pyrometry 
(TFP)  is  used  to  record  the  lempetamte  variation 
during  the  interaction  process.  A  sudden  decrease  in 
flame  temperature  intUcates  the  occurrence  local 
quenching.  The  thinning  and  disappearance  of  the  OH 
zone  during  the  quenchit^  process  is  abo  illustrated 
by  OH-PUF  imaging.  A  time-dependent  numerical 
simulation  reproduces  the  vmtex  structure  and  the 
dynamic  vortex-flame  interaction  observed  in  the 
experiment  Although  the  fast  chemistry  model  fail 
to  predkt  the  quenching,  the  calculatiao  demonstrates 
the  importance  of  tran^rt  phenomena  during  the 
vortex-flame  interaction. 

Introduction 

Lean  blowout,  high  altitude  relight  emissions,  and 
combustion  efficiency  are  of  practical  importance 
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in  the  design  of  gas  turbine  combustors.  They  all 
involve  uirbulem  combustion  {nocesses  that  need  to 
be  better  understood  if  accurate  combustor  design 
models  are  to  be  develt^ted.  Turbulent  jet  diffusion 
flames  are  used  to  gain  insight  into  turbulent 
combustion  processes  and  to  develop  numerical 
models  for  use  in  combustor  design.  Because  of 
their  i»actical  inqiottance,  they  have  been  studied  few 
over  40  years.  Although  considerable  progress  has 
been  ma^  in  understanding  turbulent  jet  diffusion 
flames,  there  still  remain  many  important  questions 
about  the  turbulent  combustion  processes  and  the 
simplifying  assumptions  retpiiied  to  effectively  model 
them.  The  coocqxual  view  of  turbulent  reaction 
zones  and  the  inieiplay  between  turbulence  and 
chemistry  are  germaiK  to  the  modeling  q)proeches. 
Recent  research  has  focused  on  whetho’  turbulent 
flames  can  be  viewed  as  an  ensemble  of  strained 
laminar  flamelets  or  distributed  reaction  zones  (Bilger, 
1989;  Correa,  1992).  In  this  paper  the  statistical 
fabric  of  a  turbulent  reaction  zone  is  viewed  as  an 
uisemUe  of  vortex-flame  interactions.  The  statistics 
of  the  problem  are  eliminated  by  isolating  and 
studying  a  single  vortex-flame  interaction  event  that 
has  similar  characteristics  to  those  occurring  in 
turbulent  jet  diflusion  flames. 

Many  of  the  basic  characteristics  of  turbulent 
jet  diffusion  flames  were  described  by  Hottel  and 
Hawthorne  (1949),  Scholefield  and  Garside  (1949), 
and  Wohl  et  al.  (1949).  According  to  their 
description,  a  vertically-mounted  jet  diflusion  flame  is 
considered  fiilly  turbulent  when  its  entire  surface  area 
becomes  a  flame  brush,  which  means  that  the  flame 
has  a  highly  wrinkled,  bumpy,  or  rough  appearance. 

Improved  visualization  techniques  have 
provided  additional  insights  into  the  structure  of  a 
turbulent  flame  and  the  events  that  happen  as  the 
flame  becomes  turbulent  (Roquemore,  et  al.,  1988). 
Although  the  turbulent  transition  processes  can  be 


observed  in  flames  established  with  fuel  tubes 
(Takeno  and  Kotani,  1977;  Takahashi.  et  al.  1982). 
they  are  more  clearly  observed  when  the  jet  flame  is 
est^lished  by  a  contoured  nozale  with  a  flat  velocity 
exit  profile  (Yule  et  al.  1981;  Eickhoff,  1982:  Savas 
and  Gollahalli.  1986;  Coats  and  Zhao  1988;  and 
Roquemore  et  al.  1988).  Al  transitional  jet 
velocities,  well  organized,  axisymmetric  vortex 
structures  are  observed  to  exist  inside  the  laminar 
flame  sur&ce  for  many  jet  duuneteis  downstream.  As 
the  jet  velocity  increases,  the  iaoer  vortices  begin  to 
coalesce  at  a  ^wnstream  location  in  the  flame.  The 
coalescence  process  is  highly  nstable  and  results  in 
the  break-down  of  the  organized  vortices  into  smaller 
but  identifiable  three  dimensiaaal  fluid  elements.  The 
flame  becomes  turbulent  when  these  small  three- 
dimensional  structures  interact  with  the  flame  creating 
die  turbulent  flame  brush.  The  hei^  at  which  this 
takes  place  is  called  the  break-point  As  the  jet 
velocity  increases,  the  break-paiat  moves  toward  the 
nozzle  exit  When  the  break-paiat  is  occurred  near 
the  Up  of  the  nozzle,  the  flame  is  considered  to  be 
fully  turbulent  Thus,  the  wrMles  or  bumps  in  the 
flasM  are  tocaUzed  prottusioas  of  the  flame  surface 
that  result  when  three-dimensianal  fluid  elements  of 
different  size,  shape,  velocity,  and  rotatkioal  strength 
interact  with  reaction  zones  of  (Hfaem  diicknesses. 

There  is  a  question  about  whether  or  not  the 
smaO  tfaree-dimensioaal  fluid  elements  Conned  in  the 
coaleaoeace  process  are  vortices,  in  a  vortex,  the  fluid 
mass  moves  around  a  common  axis  (Lugt  1979). 
Some  of  the  smaU  structures  in  turbulent  jet  flames 
are  easily  recognized  in  a  two-dimensional 
visualizations  as  vortex  pairs  that  are  associated  with 
mushroom-shaped  vortices.  Odter  fluid  dements  are 
not  so  easily  recognized  as  vortices  in  the  laboratory 
reference  frame.  However,  three-dimensional 
calculations  show  that  the  mudl  fluid  elem«?vits 
resulted  from  the  coalescence  of  two-dimensional 
vertices  are  vortex  structures  CKnio  and  Ghooiem, 
1989).  Since  this  process  is  siaiilar  to  that  occurring 
in  jet  flames,  it  will  be  assmned  that  the  three- 
dimensional  fluid  elements  that  interact  with  the 
flame  are  vortices.  The  interaction  processes  will  be 
referred  to  as  vortex-flame  interxtions. 

Individual  vortex-flame  interactions  are 
important  to  understand  for  several  reasons.  First, 
they  can  increase  or  decrease  the  turbulent  reaction 
rate.  When  the  vortex-flame  interactions  is  moderate, 
they  increase  the  surface  area  of  tie  flame  which  leads 
to  an  overall  increase  in  the  global  reaction  rate  and  a 
reduction  in  the  length  of  the  flame.  However,  if  a 
vortex  has  a  sufficiently  high  radial  velocity,  it  can 
pass  through  the  flame  creating  a  localized  hole  in 


which  no  chemical  reactions  take  place  (Takahashi 
and  SdunoU.  1990;  Chen  and  Goss,  1991;  Takahashi 
and  Goss,  1992).  The  holes  can  cause  the  flame  to 
split  or  lift  when  they  form  near  the  jet  exit.  If  the 
holes  covers  two  much  of  the  flame  surface  the  global 
reaction  rate  can  decrease.  This  can  lead  to  reduction 
in  combustum  efficiency  and  blow-ouL  Second, 
flame  vortex  interactions  may  be  building  blocks  for 
statistical  theories  of  turbulent  flames.  It  is 
ctmceivable  that  an  ensemble  of  only  a  relatively  few 
basic  types  of  vortex-flame  interactions  need  to  be 
considered  to  accurately  represent  the  statistical  result 
of  what,  on  first  sight,  mi^t  appear  to  be  an  endless 
number  of  interactions.  If  this  idea  is  supported  by 
experimental  results,  it  might  be  feasible  to  construct 
a  statistical  turbulent  combustion  model. 

Vortex-flame  interactions  have  been  studied 
experimentally  as  well  as  computationally.  Hbwever, 
most  of  the  studies  are  for  prmixed  systems  (Roberts 
et  at,  1992;  Rutland  and  Ferziger.  1989;  Poinsot  et 
aL.  1991).  Hie  studies  of  diffusioo  flames  appear  to 
involve  vortex-flame  interactions  that  ap|dy  u>  the 
flame  inieracting  with  the  large  buoyancy  induced 
vortex  stractnres  formed  outside  uansitional  jet 
flames  (Kangozian  and  Marble  1986;  Laverdant  and 
Candel  1989).  In  turbuleat  jet  diffusion  flames 
surrounded  by  nearly  still  air,  the  vortices  must  be 
transported  firm  the  tiiel  side  of  the  flame  towards  the 
diffbiian  flame  surface.  To  oiv  knowledge,  only  one 
individnal  vortex-flame  interaction  study  has  been 
performed  under  these  conditions  and  it  was 
complicated  by  the  coupling  of  a  large,  buoyancy 
indu^,  vortex  that  was  outside  of  the  flame  surface 
with  the  vortex  inside  the  flame  surface  (Lewis, 
1988). 

This  piqier  presents  the  initial  results  of  a 
vortex-flame  interaction  that  has  some  features  like 
those  observed  in  uirbulent  jet  diffusion  fUunes.  A 
counter  rotating  pair  of  ring  vortices  are  formed  by 
driving  the  flow.  Planar  visualization  and  line- 
tenqieiaiHre  measurement  techniques  are  phase-locked 
to  the  driving  pulse  so  that  the  interaction  can  be 
strobe  at  different  time  steps.  A  direct  numerical 
simulation  with  infinitely  fast  chemistry,  unity 
Lewis  number,  and  temperature  and  concentration 
dependent  transport  coefficients  (Katta  et  al.,  1992)  is 
us^  to  investigate  the  dynamic  characteristics  of  the 
vortex-flame  interaction. 

Experimental  Descriptions 

The  vortex-flame  interactions  are  produced  in 
an  externally  driven  methane  jet  diffusion  flame 
shown  schematically  in  Fig.  1.  This  is  a  very  flexible 
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system  that  allows  the  frequency,  intensity,  and  shape 
of  the  driving  pulse  10  be  varied.  A  sqwaker  connected 
to  a  25  mm  diameter  fiiei  tube  is  used  to  periodically 
generate  the  voftkes.  The  speaker  is  pla^  inside  a 
chamber  with  the  speaker  diaphragm  oriented 
vertically  upward  along  the  flow  direction.  The 
chamber  is  sealed  from  the  surrounding  air.  The  fuel 
enters  the  speaker  chandxr  through  side  inlet  ports. 
Another  fuel  inlet,  at  a  location  just  above  the 
chamber,  is  used  for  seeding  in  visualization.  A 
power  amplifier  in  conjunction  with  a  pulse  generator 
are  used  to  drive  the  speaker,  which  forces  the  fuel 
through  the  fuel  tube.  Vortices  are  formed  inside  the 
flame  near  the  tip  of  a  10  mm  diameter  long-tapered 
nozzle.  The  strength  and  the  frequency  of  the 
generated  vortices  can  be  adjusted  by  the  pulse 
generator.  The  flame  is  shielded  from  room  air 
disturbances  by  a  slowly  co-flowing  anmilus  air.  The 
air  duct  (152  mm  diameier)  consist  of  diffusion  and 
contraction  sections  and  layers  of  screens  and 
honeycombs  for  flow  straightening  so  the 
surrounding  air  flow  is  relative  quite.  The  vortex- 
flame  interaction  is  stndied  using  visualization, 
velocity  and  temperature  measurement  techniques  that 
ate  phamd-locked  to  the  driving  pulse. 

Renerive-Mie-Scanering  Imaging 

The  detailed  structures  of  the  driven  vortices 
and  the  flame  ate  visualized  using  a  Reactive-Mie- 
Scattering  (RMS)  technique  (Chen  and  Roquemote. 
1986).  The  output  beam  (532  nm)  of  a  pulsed 
Nd:YAG  laser  is  used  to  form  a  laser  sheet  (I  mm 
thickness)  vertically  along  the  center-line  of  die  jet 
flame  as  shown  in  1.  The  vortical  structures  are 
visualized  by  the  scattering  of  micron-size  Ti02 
particles  formed  from  the  chemical  reactkm  between 
the  T1CT4  vapw  seeded  in  the  fuel  and  tbecombusdcn 
product,  H2O.  A  35  mm  camera  placed  at  a  ri^t 
angle  to  the  laser  sheet  is  used  to  photograph  the 
vortical  structure  and  the  flame.  A  gated  mechanical 
shutter  mounted  on  the  camera  is  triggered  externally. 
In  order  to  visualize  the  detailed  dynamics  of  the 
flame-vortex  interaction,  a  synchronization  circuit  is 
used  to  trigger  the  Nd:Y  AG  laser,  speaker,  and  camera 
shutter.  A  master  pulse  generator  in  combination 
with  a  delay  circuit  is  used  to  create  an  adjustable 
time  delay  between  the  speaker  and  the  Nd:  Y  AG  laser 
for  phase-locked  visualization. 

Joint  RMS-PLIF  Imaging 

The  vortex  structure  and  the  OH  zone  are 
visualized  simultaneously  using  Reactive-Mie- 


Scattering  (RMS)  and  Planar-Laser- Induced- 
Fluorescence  (PLIF).  The  experimental  setup  of  the 
RMS-PLEF  system  is  shown  in  Fig.  2.  The 
frequency-doubled  laser  beam  (532  nm)  of  a  pulsed 
NtkYAG  laser  is  used  to  pump  a  dye  laser.  The 
output  UV  beam  (282.7  nm)  from  wavelength 
extender  is  uined  to  excite  the  OH-radical  (Hanson, 
1986).  A  quartz  focusing  lens  and  a  cylindrical  lens 
are  combirt^  to  expand  the  UV  laser  sheet  vertically 
along  the  center-line  of  the  flame.  A  gated 
intensified-CCD  camera  is  used  to  record  the 
fluorescent  ima^  of  the  OH-radical  at  a  right  vigle  to 
the  truer  sheet.  The  phase-locked  measurement  is 
accomplished  by  precise  synchronization  between  the 
laser,  speaker  pulse,  and  CCD  camera.  The 
intensifier,  gated  by  a  pulse  generator  with  1.5  ms 
gate  width,  is  used  to  diminate  the  flame  background 
emission.  The  hfie  scattering  of  the  T1O2  panicle 
marks  the  r^ioii  where  the  fuel  and  products,  H2O, 
are  mixed. 

Uiicr-Pnwiicr-AnciMii^ 

A  frequency-shifted  Laser-Doppler- 
Anemometry  (LDA)  system  is  used  for  velocity 
measurement  The  fuel  is  seeded  with  AI2O3 
particles  (5  pm  mean  diameier).  The  time  history  of 
the  velocity  is  measured  by  employing  a  phase- 
locked.  conditumal  sanqrling  method.  A  5  ms  time 
window  is  used  to  register  KXX)  samples  at  each 
phase  angle. 

Thin-Filament  Pviometrv 

The  line  temperature  is  measured  by  Thin- 
Filament  Pyromeiry  (TFP)  technique  (Vilinqxx;  and 
Goss.  1988)  as  shown  schematically  in  Fig.  3.  A 
fine  SiC  fiber  with  approximately  15  pm  diameter  is 
supported  and  placed  horizontally  across  the  center  of 
the  flame.  The  thermal  radiation  from  the  heated 
filament  is  focused  onto  a  ten-facets  spinning  mirror. 
An  InGaAs  detector  with  spectral  range  of  900  to 
1500  nm  is  used  to  measure  the  radiative  signal  from 
the  filament  The  emission  intensity  along  the 
filament  is  scamied  continuously  onto  tlw  detectin'  by 
the  ginning  minor.  A  12-bit  waveform  recorder  is 
used  to  digitize  and  store  fllament  intensity.  AHe-Ne 
laser  and  a  photodiode  are  used  as  the  triggering  signal 
for  initiating  data  acquisition.  The  line  intensity  of 
the  fllament  is  scanned  on  to  the  InGaAs  detector  at  a 
2  kHz  frequency.  The  InGaAs  detector  signal  is 
digitized  by  the  waveform  recorder  at  5  MHz 
frequency.  A  program  is  used  to  interface  the 
waveform  recorder  and  a  Macintosh  Ilfx  for  data 


3 


acquisitk)ii  and  tmsfer.  The  digitized  1TP  dtta 
consist  of  1024  points  in  spatial  diiectian  and  S12 
segments  in  temporal  direction.  The  reconfed  TFP 
data  has  spatial  and  temporal  resoliniaos  of  33  lU" 
and  300  |is,  respectively.  The  time  trace  of  the 
driving  pulses  is  digitized  by  a  muld-channcl  16-bit 
A/D  converter.  The  TFP  data  is  collected  at  20 
different  axitd  locations  above  the  nozzle  exit  The 
flame  background  emission  is  also  recorded  at  those 
same  locations.  Background  subtraction  with  precise 
phase-matching  is  performed  or  the  TFP  to 
eiiminalB  the  flame  emission.  The  recorded  time  trace 
of  the  driving  pulses  ensure  precise  phase-msiching  in 
datt  reduction.  Temperature  is  cakiilaind  from  the 
processed  TFP  data.  The  detailed  calibraiion  sod 
correction  of  the  TFP  measurement  has  been 
discussed  by  Vilimpoc  and  Goss  (198S).  Phase- 
avenging  on  temperature  is  performed  over  7  cycles 
at  each  axial  location.  The  phase-averaged 
teoqierature  profiles  calculated  at  each  location  are 
used  to  recoMtnict  the  two-dimeosiorial  temperamre 
field  m  different  phase  angles. 

Nupyw*!  Srfteme 

The  difibsian  flame  coosideted  in  this  work  is 
formed  between  a  central  methane  jet  and  a  low  speed 
co-annolar  air  flow.  Time  dependent  governing 
equations,  expressed  in  cylindiicalcoor«Biiates.consiat 
of  mass  continoity,  axial  and  radial  moatentnm 
conservation,  and  two  scalar  conservation  equations 
(Katta  et  al.  1992).  Shvab-ZeTdovich  formntatiaa  in 
conjunction  with  the  flame-sheet  modd  is  assumed, 
This  formulation  is  primarily  used  in  undentantBng 
the  mechanism  leading  to  large  vortices  that  are 
observed  in  the  experiments  and  their  effects  on  the 
flame  surface.  Body-force  term,  doe  to  the 
gravitational  field,  is  included  in  the  axial  momemnm 
equatitm.  The  following  global  chemical  kineiics 
between  methane  and  oxygen  is  coosideted : 

CH4  +  2  02  +  7.3  N2  ->  2  H2O  +  CO2  +  7.3  Ni 

Nitrogen  in  this  system  is  assumed  to  be  an 
inert  gas.  The  system  of  governing  equations  is 
completed  by  using  the  state  equation.  Tcansport 
properties  (viscosity  and  diffusion  coefficient)  are 
considered  to  vary  with  temperature  and  species 
concentrations.  The  enthalpy  of  each  species  is 
calculated  from  4>olynomial  curve-fits,  wlule  the 
viscosity  of  the  individual  species  is  estimaied  from 
Chapman-Enskog  collision  theory.  The  binary 
diffiision  coefficient  between  any  two  ^lecies  on  the 
fuel  side  of  the  flame  is  assumed  to  be  identically 


equal  to  that  of  the  fuel  and  nitrogen.  Similarly,  on 
the  oxidizer  side  of  the  flame  it  is  made  identical  to 
that  of  the  oxygen  and  nitrogen.  Chapman-Enskog 
theory  and  the  Lennard-Jones  potentials  have  been 
used  to  estimate  these  two  binary  diffusion 
coefficients. 

The  finite-difference  form  of  the  governing 
equations  are  constructed  on  a  staggered  grid  system 
baaed  on  an  implicit  QUICKEST  numerical  scheme. 
It  is  third-order  accurate  in  both  space  and  time  and  has 
a  very  low  numerical  diffusion  error.  At  every  dme- 
stq),  the  pressure  field  is  accurately  calculated  by 
solving  the  system  of  algebraic  pressure  Poisson 
equatioos  simultaneously,  at  all  grid  points  by  using 
LU  decomposition  technique.  An  orthogonal  grid 
system  with  rapidly  expanding  cell  sizes  in  both  z  and 
r  directions  is  utilized.  The  outer  boundaries  of  the 
computational  domain  are  shifted  sufficiently  far 
enoi:^  to  minimize  the  propagation  of  disturbances 
into  die  region  of  interest  Outflow  boundary  in  these 
flows  is  the  most  difficult  one  to  treat  berause  the 
flow  leaving  this  boundary  continuously  evolves  in 
time  as  the  large  vortices  convect  trough  this 
boundary.  A  sim|de  extrapolation  procedure  with 
weighted  zero- and  first-older  terms  is  used  to  estimate 
the  flow  vari,‘ibles  on  the  boundary  in  the  present 
caicniatinos.  Hie  main  criterion  used  in  selecting  the 
weighting  functions  is  that  the  vortices  crossing  the 
outflow  boundary  should  leave  smoothly  without 
being  (Sstorted. 

Calculations,  using  a  mesh  size  of  ISl  x  61 
have  been  made  for  the  unsteady  flow  in  a  230  mm  x 
ISO  mm  computational  domain.  Flat  initial  velocity 
pnrfiles,  with  the  measured  mass-averaged  values  of 
U  and  0.4  m/s,  have  been  used  at  die  exit  planes  of  a 
10  mm-cemer-nozzle  and  a  ISO  nun  annulus-duct, 
reflectively.  Because  of  the  gravity  imm  in  the  axial- 
momentum  equation  die  computed  flame  developed  a 
low-frequency  oscillation  of  about  20  Hz.  However, 
no  roil-op  (rf’  the  buoyancy-induced  vortex  occurred 
within  an  axial  distance  of  230  mm. 

The  large  driven-vortices  in  the  experiment  are 
then  simulated  by  periodically  perturbing  the  fuel 
mass-flow-rate  at  the  exit  of  the  fuel-nozzle  at  a 
frtqiroiicy  of  30  Hz.  The  modified  velocity  at  this 
section,  for  a  period  of  one  perturbation  cycle,  is 
shown  in  Fig.  4  along  with  that  of  the  measured  one. 
Note  that  there  is  a  suction  pulse  of  -1.3  m/s  for  a 
period  of  about  6  ms  making  the  modified  velocity 
zero.  This  stagnant  condition  at  the  nozzle  exit  pulls 
the  flame  surface  and  thereby,  the  viscous  fluid  closer 
to  the  axis.  The  following  sudden  injection  of 
additional  fuel  at  a  velocity  of  6  m/s  (modified 
velocity  is  7.3  m/s)  makes  a  collision  with  the  highly 
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viscous  combustkm  pnxiucts.  This  piocess  results  in 
sending  an  axisymineiric,  counter  rotating  pair  of  ring 
voitioes  iadii%  toward  the  flame  surface.  Thegrowth 
and  the  intencdoo  of  this  vortex  with  the  flame  is 
further  discussed  in  the  results  and  discusskn  sectiaiis 

The  purpose  of  this  study  is  to  investigate 
vortex'flame  interactions  that  result  in  flame  stretch 
and  extinction  processes  thm  are  similar  to  those 
observed  in  turbulent  jet  diffusion  flames.  With  this 
purpose  in  mind,  but  realizing  that  actual  voitex- 
flaim  interactions  in  turbulent  combustion  processes 
are  three-dimensional  and  that  this  seriously 
complicates  both  the  modeling  and  experimental 
efforts,  the  following  characteristics  of  the  vortex- 
flame  interactioo  were  sought  1)  The  outward  radial 
velocity  of  the  vortex  should  be  controllable  so  that 
vorteX'flame  interactions  results  in  either  the 
production  of  a  flame  bulge,  which  does  not 
extinguish  the  flame,  or  a  flame  bulge  that  is 
stretched  to  the  pram  that  the  flame  is  extingniriied.  a 
bole  is  produced.  Similar  effects  ate  known  tooccur 
in  turfauleat  jet  diffusion  flames.  2)  The  interaction 
should  be  axisyrrunetric,  since  this  would  fipciiitawt 
the  use  of  direct  numerical  simulation  and  also 
simplify  the  experimenL  3)  The  vortex  should 
interact  with  a  tto,  laminar  flame  so  that  a  flame- 
sheet  model  can  be  used  in  the  simulation.  4)  The 
shape  of  the  vortex  was  not  a  major  considetmioo 
since  visualizatkms  have  shown  that  many  different 
shape  fluid  elements  interact  with  the  flame. 
However,  a  counter  rotating  vortex  pair  ate  oAra  seen 
in  two-dimensional  visnalizroioos  of  torbdlem  jets 
and  turbulent  jet  flames  so  it  would  be  to 

produce  this  type  of  vortex  sotictote.  S)  The  vortex- 
flame  interaction  should  be  periodically  ptoduoed  so 
that  phased-locked  measurements  could  be  auke. 
This  allows  the  vortex-flame  interaction  to  be  smdied 
at  discrete  time-delays  diving  the  interactian  process. 

A  significant  investigation  was  required  to 
detomine  the  driving  conditions  and  the  fuel  flow  rate 
that  would  generate  a  vortex-flame  interaction  diat  had 
the  characteristics  described  above.  RMS 
visualizations  w«e  used  to  systematically  tnve«igate 
the  influences  of  frequency,  intensity,  and  pulse  shape 
on  the  vortex  formation,  growth,  and  flame 
interaction  processes  for  different  air  and  fuel  flow 
rates.  Several  problems  were  encountered.  Hrst,  it 
was  difficult  to  find  conditions  at  which  a  vertex 
would  grow  radially  to  the  extent  needed  to  imeract 
with  the  flame  surface.  Second,  when  such 
conditions  were  found,  the  flame  would  usualiy  lift. 


hence  (he  interaction  process  would  not  be  periodic. 
Third,  the  flame  would  go  three-dimensional.  Fourth, 
large  buoyancy  induced  vortices,  formed  outside  the 
flame,  would  couple  with  the  inner  vortices  formed 
by  the  driving  pulse,  so  the  system  would  become 
very  compUcated.  It  was  found  that  large  internal 
vortex  structures  were  most  evident  for  flow  rates 
below  about  2  m/s  and  at  driving  frequencies  of  30. 
4S.  60.  and  90  Hz.  The  calculated  resonance 
frequency  of  the  fuel  tube  (or  the  organ  pipe  effect) 
was  240  Hz.  The  fretpiencies  of  rapid  vortex  growth 
were  subharmonics  of  this  frequency,  or  linear 
combinations  of  the  subharmonics.  However, 
increasing  the  amplitude  of  a  sinusoidal  wave  at  a 
favorable  fiequency  did  not  produce  a  desirable  vortex- 
flame  interaction  for  the  reasons  stated  above.  The 
key  parameta  that  lead  to  the  desired  vortex-flame 
inteiactkn  was  the  shape  of  the  driving  pulse. 

A  suiiable  vortex-flame  interaction  is  produced 
when  the  fud  flow  rate  is  maiiitain<»rf  at  a  mean  nuxs 
velocity  of  1.S9  m/a  and  an  annulus  air  velocity  of 
0.39  m/s.  Without  driving,  the  flame  is  laminar  up 
to  a  height  of  SO  mm  above  the  nozzle  dp.  When  a 
fast  rising  pulse,  with  a  12  ms  exponential  decay,  is 
applied  to  the  qpeaker  at  a  rate  of  30  Hz.  a  counter¬ 
rotating  vortex  ting  is  obaenred  to  produce  a  balge  in 
the  Uue  flame  surfoce.  Increasing  the  iniensior  of  the 
driving  pulse  causes  a  viable  hole  to  form  in  the 
stretched  flame  at  a  height  of  20  mm  above  the 
nozzle.  The  vortex-flame  interaction  is 
axisymmetric.  periodic,  md  very  rqteatable.  The 
dynamics  of  vortex-flame  interaction  can  be  observed 
with  0.5  ms  resolution  of  die  time  delay. 

The  sharp  rise,  12  ms  decay  driving  pabe  to 
the  speaker  results  in  the  axial  vehxify  pulse  Aown 
in  F^.  4.  The  Phase-locked  velocity  measurements 
were  made  on  the  nozzle  centerline,  at  z  «  4  mm.  A 
thousand  velocity  realizations  were  averaged  at  each  1 
mstimedelay.  Data  were  collected  at  33  time-delays 
which  coneqxmds  to  one  period  of  the  30  Hz  driving 
frequency.  The  time-delay  was  measured  from  the 
leading  edge  of  the  driving  pulse  to  the  middle  of  the 
LDA  data  sampling  window.  The  sampling  window 
had  to  be  opened  to  S  ms  to  obtain  a  reasonable 
sampling  rate.  This  means  that  the  measurements  at 
each  time-delay  were  actually  avmaged  over  two  and  a 
half  time-delays  on  each  side  of  the  recorded  time- 
delay.  This  resulted  in  the  large  uncertainties  noted 
by  the  error  bars  on  the  velocity  data. 

It  is  seen  from  Fig.  4  that  the  mean  axial 
velocity  increased  from  U  m/s  and  reached  a  peak  of 
7  m/s  in  about  8  ms.  This  was  followed  by  a  gradual 
decay  to  1.5  m/s  over  a  period  of  22  ms.  There 
appears  to  be  three  inflection  points  in  the  dec^  pan 
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the  velocity  pulse.  They  are  locaiBd  at  t «  16, 22 
and  28  ms.  Thw  atfleciian  poimi  me  probably  due 
to  the  mechanical  dmping  of  the  speaker. 

Phase-Locked  RMS  visualizations  at  8.  10. 
12  and  14  ms  time-delays  are  shoam  Cnm  left  to  right 
in  Fig.  S(a).  The  vortex  structures  cqNured  by  the 
RMS  visualizations  can  be  recognized  by  their  dark 
green-yellow  appeatmice.  These  irnagBS  are  ftoaen  by 
the  10  ns  flash  the  laser  sheet  passaig  through  the 
center  of  the  flmne  and  thus,  aocmaiely  repreaem  the 
time  mid  spatial  evolution  of  the  vortices.  This  is 
not  necessarily  the  case  for  the  images  of  the  flame. 

The  apparent  flame  locations  in  Fig.  S(a)  are 
marked  by  the  yellow  and  blue  colon  in  the 
photograph.  The  yellow  color  is  lirom  the  blacfcbody 
radiation  of  the  soot  particles  and  the  bine  is  the 
emissions  from  the  CH  radicals  in  reaction  zone  of 
the  flame.  Both  the  yellow  and  Mae  flames  are  three- 
dimensional  structwes  that  are  sapetiaqnaed  on  to  a 
two-dimensional  phoiogn#.  However,  m  dm  outer 
radial  location  of  the  flame,  they  can  be  a  reasonable 
two-dimensional  representations  of  the  flame  surface. 
Unfortunamly  the  camera  shutter  had  10  be  opened  for 
S  ms  to  capture  enough  light  to  see  the  Mae  flame. 
This  means  that  the  flame  images  me  not  truly 
leilecsing  the  flmne  *10001100"  at  the  icpoited  time- 
delay.  This  is  evident  from  the  hnage  of  dm  Mue 
flame  in  the  8  ms  time-delay  cam.  For  aumple, 
there  is  an  appaiem  hole  in  the  Mae  flaom  that  is 
located  above  the  vertex.  This  is  a  frrise  oaage  and 
only  reflects  the  fact  that,  during  the  S  ms  shmter 
opening,  centered  about  the  8  ms  time-delay,  the  Mue 
flame  spent  most  of  the  time  at  the  locatioa  where 
the  bde  was  formed.  However,  m  dm  higher  time- 
delays,  the  location  of  the  blue  flmne  doesn't  appear 
to  change  much,  which  indicaies  that  the  actual 
changes  in  the  location  of  the  bhm  flame  during  the 
shutter  opening  are  reasonaMy  conect 

The  vortex-flame  inieractian  are  evident  in  Fig. 
S(a).  At  the  8  ms  time-delay,  a  huge  vortex  that 
rotates  downward  is  observed.  There  is  also  a  counter 
rotating  "secondary”  vortex  on  tq>  of  the  large 
"primary"  vortex  but  it  is  to  snudi  to  be  clearly 
observed.  Two  milliseconds  later,  however,  the 
"secondmy"  vortex  is  as  large  as  the  "piimaty*  vortex. 
It  is  clear  that  they  form  a  pair  of  counter  rotating 
vortex  rings.  Almost  all  of  die  oittwaid  radial  motion 
of  the  vortices  occurs  between  0  and  10  ms.  This 
corresponds  to  the  time  that  the  axial  vdodiy  teaches 
its  peak,  as  notett  in  Fig.4.  The  voitex  interaction 
with  the  flame  produces  a  hole  that  occurs  between  8 
and  10  ms.  This  is  confirmed  by  the  OH  images  in 
Fig.  9.  It  is  also  confirmed  in  Fig.  9  that  the  vortex 
never  extends  radially  beyond  the  flame  surface.  This 


indicates  dim  flame  stretch  is  involved  in  the  flame 
quenching  process.  Once  the  hole  is  formed  and 
vortex  has  reached  its  maximum  radial  location.  The 
vortex  pair  then  appear  to  rotate  inward  as  can  be 
observed  for  time-delays  of  12  and  14  ms.  The  hole, 
which  is  clearly  a  ring,  as  noted  by  the  dark  band 
across  the  flame,  remains  open  for  a  long  time. 

The  relative  axial  velocities  of  the  fuel  at 
different  radial  locations  are  reflected  by  the  line 
connecting  the  vortices  on  the  tight  and  left  sides  of 
the  flame  in  Fig.  S(a).  The  center  of  the  connecting 
line  probably  represenu  a  stagnation  point  as 
observed  by  Lewis  et  al.  (1988).  At  the  8  ms  time- 
delay.  the  vortex  connecting  line  is  curved  upward, 
indicating  that  the  vortices,  at  their  outer  radial 
location,  have  a  higher  axial  velocity  than  the  fluid  at 
the  center  of  the  jet  At  10  ms,  the  axial  vekicities, 
at  die  height  of  the  vortices,  are  nearly  independent 
of  radial  location.  For  the  12  and  14  ms  time-delays, 
the  central  region  of  the  flow  is  moving  much  faster 
than  the  outer  radial  region  where  the  vortices  are 
located.  Indeed,  the  stagnation  point  moves  upward 
about  10  mm  between  10  and  12  ms  time  ^lays. 
This  give  a  velocity  of  S  m/s,  which  is  in  agreement 
with  the  data  in  Hg.  4  for  the  1 1  ms  time-delay. 
There  is  vesy  Bitle  vertical  motion  of  the  vortices  at 
thne-dehqrs  between  10  and  14  ms;  however,  there  is 
a  very  large  aocekrarian  of  the  central  jet  fluid  above 
tbevortioes.  This  is  caused  by  inward  radial  motion 
of  the  vottioes  which  produces  a  jetting  action  of  the 
central  fluid.  It  is  anticipated  that  this  will  suck  air 
in  towards  the  fuel  stream  to  replace  the  volume  of 
fluid  being  aoederated  upwards.  The  accelerating  fluid 
cannot  be  reidaced  by  tte  fluid  below  the  vortex  pair 
because  its  axial  velocity  is  decreasing  (See  Fig.  4  for 
time-delays  greater  than  10  ms)  by  the  inward 
contraction  of  the  vortices. 

The  numerical  simulation  at  approximately  the 
same  thne-ddays  as  those  of  the  experiments  of  Fig. 
S(a)  are  shown  in  Fig.  S(b).  Composite  figures  of 
isotherms  (left)  and  color-coded  particle  trackings 
(right)  are  presented  for  each  time-dday.  Also  shown 
in  the  flgure  are  the  predicted  instantaneous 
stoichiometric  mixture  fraction  locations,  marked  by 
the  white  dots  outside  the  particle  trackings.  The 
numerical  simulation  reproduces  the  evolution  of  the 
primary  vortex  of  the  experiments  using  the  particle 
tracking  technique.  The  counter-rotating  vortex  ring 
is  also  evident  in  the  simulation.  The  isotherm 
representation  did  not  identify  the  vortex  pair  that  is 
shown  by  particle  trackings.  In  fact,  the  temperature 
field  prediedon  seems  to  suggest  the  existence  of  a 
single  vortex  region.  The  numerical  simulation 
appears  to  produce  nearly  the  same  time  and  spatial 


evolution  of  the  ring  vortex  as  observed  in  with  the 
experimental  visualizations. 

The  simulation  shows  a  thinning  of  the 
reaction  zone  aid  a  decrease  in  the  flame  tempeiaiure 
as  the  flame  is  being  stretched  by  the  vortex-flame 
interaction.  This  is  noted  by  the  temperature 
contours  in  Fig.  S(b).  The  thinning  of  the  flame 
seems  reasonable  but,  the  actual  reduction  in  the 
stoichiometric  flame  temperature  is  not  physically 
correct  for  a  flame  sheet  model  with  fast  chemistry. 
This  happens  because  the  reaction  zone  becomes 
much  thinner  than  the  grid  spacing.  The  grids 
containing  the  stoichiometric  reaction  zone  also 
contains  an  excess  of  fuel,  so  the  inadequate  grid 
resolutioa  results  in  a  lower  flame  temperature. 

Phase-locked  TFP  measurements  at  10  and  14 
ms  are  shown  in  Fig.  6.  The  TFP  temperature  data 
are  superimposed  on  their  respective  RMS  images. 
The  results  show  that  the  local  peak  tenqietaiiires  (at 
fixed  axial  heights)  are  usually  situated  along  the 
luminous  flame,  as  one  would  expect  for  a  laminar  or 
near-fauninar  jet  diffusion  flame,  except  for  the  r^ion 
where  the  vortex  is  located.  The  phase-averaged 
temperature  profiles  show  doable  peaks  in  the 
measurements  of  the  second  and  third  filaments  above 
the  burner  tube  at  10  ms  time-delay,  and  in  the  upper 
filaments  at  14  ms  time  delay.  At  14  ms  thnedelay, 
the  outer  peak  temperature  at  z  >  18  mm  (fifth 
filament  above  the  burner)  is  850  K  at  a  riulial 
location  r  a  12  mm.  This  ra^kxatioa  is  where  the 
local  quenching  of  the  "blue"  flame  is  visually 
observed. 

The  time  trace  of  the  filament  temperature  in 
Hg.  7,  illustrates  that  the  outward  movement  of  the 
flame,  as  noted  by  the  location  of  the  peak 
temperature,  is  accompanied  by  a  decrease  in  the  L  jt 
zone  thickness  and  a  reduction  in  the  peak 
temperature.  The  peak  lemperatute  location  is 
superimposed  on  the  TFP  time  trace.  When  an  flame 
bulge  is  observed,  the  "flame”  lenqMtature  decreases 
to  below  the  detection  limit  of  800  K.  This  decrease 
in  "flame”  temperature  is  accompanied  by  the 
appearance  of  a  second  tenqieratore  peak  at  a  radial 
location  closer  to  the  jet  center  line.  This  second 
peak  corresponds  to  the  inner  peak  shown  by  the  TFP 
data  in  Fig.  7.  The  radial  locations  of  the  flame  (the 
outer  peak)  and  the  inner  peak  me  also  shown  in  Fig. 
7  as  a  function  of  time  delay.  As  the  flame  moves 
outward,  the  inner  peak  in  Fig.  7  starts  to  appear. 
When  the  flame  is  quenched  by  the  vortex-flame 
interaction,  the  Inner  peak  temperature  becomes 
higher  than  the  outer  p^  temperature.  The  time- 
evolution  of  the  outer  a^  itmer  peak  temperatures  are 
shown  in  Fig.  8.  The  duplicate  sets  of  triangles 


Ulustraie  the  symmetry  of  the  TIT  measurements. 
The  outer  peak  temperature  is  seen  to  decrease  from 
1800  K  to  800  K  within  the  first  3  ms  during  the 
apUi  growth  of  the  flame  radius.  The  local  peak 
tempenttte  gradually  recovers  the  unperturbed  flame 
lonpeiM 

time  response  of  the  filaments  is 
significant  in  comparison  with  the  temperature  time 
traces  of  the  flame.  The  time  response  of  the 
filament,  as  the  temperature  is  instanuneously 
drqjped  from  1750  K  to  300  K,  is  shown  in  Fig.  8. 
The  flow  remains  at  300  K  fmr  cooling  times  of  1.  3. 
and  5  ms.  The  response  m  the  filament  for  these 
cooling  limes  and  a  gas  flow  velocity  of  1  m/s  are 
also  shown  in  Fig.  8.  The  calculated  filament 
tempetaimes  decrease  from  1800  K  to  1360  K,  850 
K.  650  K  for  1ms,  3  ms,  and  5  ms  cooling  pulse 
widths,  respectively.  A  3.5  ms  cooling  period  is 
requiied.  to  reach  the  detector  cut-off  temperature  of 
800  K.  This  suggest  that  the  filament  is  exposed  to  a 
low  tempenture  rnivironnmitt  as  shown  by  the  time 
trace  in  Hg.  7.  The  time  response  of  the  filament  is 
cunently  being  investigated  for  different  convective 
flow  condMons  (Post  and  Gcrd.  1993). 

Foot  OH-PLIF  images  at  t «  9, 11. 13.  and  15 
ms  are  Aown  in  Fig.  9.  These  LIF  images  were 
taken  with  the  intensifier  gate  width  set  at  1.5  ms. 
Also  shown  on  the  LIF  images  are  the  RMS 
visnalizmion  taken  simultaneously.  The  vortex 
iniemctian  with  the  OH  zone  (or  hot  zone)  is  clearly 
iUustraied  in  Fig.  9.  As  the  "primary”  ring  vortex 
evolves  and  grows  in  size,  the  OH  becomes 
contoured.  A  local  thirming  of  the  OH  zone  is 
evident  in  the  visualization  taken  at  the  9  ms  time- 
delay.  The  radial  movement  (or  growth)  of  the  vortex 
is  suCBdent  to  form  a  flame  bulge.  Further  outward 
radial  movement  of  the  ring  vortex  results  in  a 
decrease  in  the  local  (XI  fluorescence  signal  below  the 
detection  IkniL  At  his  point  a  hole  is  formed  in  the 
flame,  as  noted  in  the  11  ms  time-delay  image.  This 
local  "quenching"  persists  for  a  few  millisecond,  as 
illustrated  by  the  images  at  t  =  13  and  15  ms.  A 
protruding  of  the  ring  vortex  is  also  evident  at  these 
two  time  delays. 

Piscussion 

The  vortex-flame  interaction  of  Figs.  4-9 
shows  local  extinction  due  to  the  flame  stretch  as 
opposed  to  other  types  of  the  vortex-flame 
interaction,  for  example,  that  resulting  from 
penetration  of  the  fast  moving  vortices  through  the 
flame  (reaction)  zone  as  in  the  formation  of  the 
flameiets  in  jet  diffusion  flames  (e.g.,  see  Roquemore 


et  aL,  1^).  The  tenglh  acak  of  the  ring-vodex  pair 
of  the  present  experiments  is  mudi  larger  than  that  of 
tlw  flmne  zone,  and  the  tone  scak  much  longer. 
These  charactmistics.  however,  need  not  always  to  be 
the  case  for  turbulent  flaaaes  in  propulsiaa  systems. 

The  stretching  of  the  flame  due  to  the  vortex 
action  appears  to  lead  to  local  extinction  of  the  hot 
(orOH)zone.  To  obtain  a  belter  understanding  of  the 
transport  process  that  may  be  critical  to  the  vortex- 
flame  interaction  studied,  the  numerical  simulation  at 
a  phase  angk  similar  to  the  10  ms  time  delay  of  the 
experiment  is  studied.  The  results  ate  shown  in  Fig. 
10(a).  using  a  similar  composite  image  format  of 
Fig.  S(b).  The  radial  prafUesof  temperature  and  mole 
fractions  at  three  axial  kcaiioos.  z  «  18. 28  md  40 
nun.  are  shown  in  Fig.  10  (b).  by  the  side  of  the 
composite  image.  The  phaae  angle  of  10  ms  time 
delay  is  chosen  because  it  is  in  the  initial  stage  of  the 
local  extinction  and  the  coonier-rataiiog  ring  vortex  is 
just  being  formed.  Thus,  a  better  understanding  can 
be  obtained  for  the  transport  processes  critical  to  these 
phetKxnena. 

On  the  focnudon  of  the  counter-rotating  ring 
vortex,  the  stmulation  riiows  a  high  radial  velocity 
component  in  the  primary  ring  vortex.  The  particle 
trackings  clearly  show  an  outward  motion  of  the  jet 
fluid.  This  outward  motion  was  deflected  at  the  flame 
location,  due  to  its  inability  to  penetrate  the  flame 
zone,  perhaps,  due  to  the  btnkr  the  high  viscosity 
at  the  flame  andAv  the  hsge  deosi^  talk)  between  the 
entrained  fluid  and  the  flame.  Whik  the  downward 
deflected  fluid  is  being  emrained  into  the  primary  ring 
vortex,  the  iqrward  moving  fluid  forms  a  counto*- 
rotating  ring  vortex  on  top  of  the  primary  vortex. 
The  partick  tracking  of  the  radial  motion  mimics  a 
"stagnation-like"  layer  shown  by  RMS  visualization. 

On  die  transport  of  !g)ecies  and  thermal  energy, 
the  radial  proflles  of  temperature  and  species 
concentration  at  z  «  28  mm.  Le..  location  inside  the 
vortex,  show  sharp  gradients  of  these  quantities  by 
both  sides  of  the  stokhionietiic  location.  These 
sharp  gradients  suggest  an  enhanced  transport  of 
reactants  (as  can  be  measured  by  local  scalar 
diashwtion  laies)  into  the  reaction  zone,  and  enhanced 
conductive  heat  loss  from  the  flame  zone.  When  the 
fast  chemistry  assumption  is  employed,  the  prediction 
stales  that  all  the  reactants  tranq;ioned  to  the  reaction 
zone  will  be  consumed  by  the  chemical  reaction,  and 
at  the  unity  Lewis  number  limit,  the  increased  heat 
release  rue  due  to  enhanced  diemical  reaction  rates 
will  be  balanced  by  the  enhanced  conductive  heat  loss 
as  reflected  by  thc'sharp  temperature  gradients.  Thus, 
an  adiabatic  flame  is  always  present  at  the  predicted 
stoichiometric  location.  In  reality,  however,  the 


reaction  rate  can  not  keep  up,  without  limit,  with 
enhanced  transport  of  the  leacianis.  As  a  result,  local 
extinction  will  occur.  In  fact,  the  phase-locked  data 
taken  at  10  ms  time  delay  (through  IS  ms  time  delay) 
suggests  that  local  extiiKtion  indeed  is  present  An 
interesting  observation  of  the  predicted  temperature 
profiles  is  that  at  z  s  28  mm  of  Fig.  10(b),  the 
species  concentration  gradient  is  so  sharp  that  the 
present  numerical  grid  resolution  is  not  sufrkient  to 
capture  the  instantaneous  stoichiometric  location. 

On  the  instantaneous  temperature  proflles. 
doubk  peaks  in  temperature  are  predicted  at  the  axial 
location  where  the  primary  vortex  is  present,  e.g..  z  = 
28  mm.  The  prediction  shows  that  the  outer  peak  at 
this  location  (z  =  28  mm)  corresponds  to  the 
stoichiometric  location;  thus,  it  has  a  higher 
temperature  than  the  inner  peak.  The  TFP  data  of 
Rgs.  6(a).  however,  show  that  the  maximum 
temperature  appears  at  the  inner  peak  location  rather 
than  the  outer  peak.  The  pr^cted  inner  peak 
temperature  m  z  »  28  mm  is  a  result  of  the  convective 
transport  of  hot  combustion  products  into  the  primary 
vortex.  At  the  fast  chemistry  limit,  the  predicted 
inner  peak  temperature  is  S(X)K.  The  TFP  data  show 
an  inner  peak  temperature  of  1300  K  and  the  outer 
peril  temperature  below  ^  K.  The  convective 
tmqxMt  of  combustion  {uoducts  from  an  iqistream 
flame  bcaiion  can  not  explain  the  much  high^  inner 
peak  temperature  determined  by  the  TFP 
measurements.  Due  to  locri  quenching,  it  is  likely 
that  oxygen  can  be  entrained  imo  the  vortex  core  or  it 
can  leA  through  the  reaction  zone.  As  a  result, 
partially  premixed  flames  can  exist  inside  the  ring 
vortex,  as  well  as  the  convoluted  stokhiomeuic 
surfaces  To  verify  the  aforementioned  hypotheses  for 
the  temperature  dwUe  peaks,  the  experiments  call  for 
concentration  measurements  inside  the  ring  vortex, 
and  to  advance  our  understanding  on  the  type  of  the 
vortex-flame  interaction  studied,  the  numerical 
simulation  requires  to  incoipOTate  the  finite  rate 
chemistry  for  study  of  the  reaction  zone  structures 
inside  the  vortex  and  the  extinction  phenomenon. 

Summary  and  Conclusion 

This  paper  reports  the  initial  results  of  a  joint 
expetimenud  and  numerical  study  of  the  vortex-flame 
interaction  in  difliision  flame. 

The  experimental  system  periodically 
reproduces  an  isolated  vortex-flame  interaction  event 
that  [xovides  a  unique  opptMtunity  for  study  of  a 
specific  type  of  the  vortex-flame  interaction,  that 
leads  to  local  quenching  of  the  flame.  The  time 
evolution  and  the  dynamics  of  the  vortex-flame 


8 


% 


interaction  are  quantified  by  phase- locked 
measuiemenis  of  the  ftiel-jet  exit  velocity  and  the 
now  field  tempeiatuie.  along  adth  |diaae-locked  RMS 
visualization  and  PUF  OH  imaging. 

A  ring  vortex  is  periodic^y  reproduced  by  the 
drivmg  of  the  fuel  jet  using  a  acoustic  speaker, 
starting  near  the  fuel-jet  exit  location.  This  vortex 
grows  in  size  as  time  evolves  and  pushes  the  location 
of  the  luminous  flame  radially  outward.  A  flame 
bulge  is  subsequently  formed.  During  the  growth  of 
the  primary  vortex,  e.g.,  vortex  entninmem  of  the 
fuel  jet  fluid,  and  the  time-varying  mass  addition  of 
the  driving,  a  counter-rotating  vortex  ring  (secondary 
vortex)  is  formed.  This  secondary  vortex  was 
initially  situated  on  top  of  the  primary  vortex.  A 
"solid-body'*  type  of  rotation  is  subsequently  observed 
from  RMS  visualizations,  showing  the  rotation  of 
the  secondmy  vtxtex  ring  around  the  primary  vortex, 
or  vice  versa.  This  "solid-body"  roution  is  a  lesult  of 
the  faster  moving  primary  vortex. 

When  the  flame  bulge  is  formed,  the  PLIF 
images  shows  a  local  thinning  of  the  OH  zone  at  the 
maximum  radius  location.  The  thinning  of  the  OH 
zone  suggests  the  flame  being  stretched  due  to  the 
interaction  of  the  vortex  and  the  flame.  The 
conliiHMd  growth  of  the  vortex  size,  or  the  flame 
bulge  radius,  eventually  results  in  a  local 
disappearance  of  the  OH  flumescence.  coincideat  with 
whm  the  filament  emissive  power  falls  bdow  the 
detection  limit  (i.e.,  800K).  At  this  comfition  the 
TFP  data  showed  double  peaks  in  radial  profiles, 
onlfte  without  artificial  driving,  the  TFP  data  show  a 
single  peak  located  outside  the  luminous  flame.  The 
double  peak  temperature  profile  shows  a  higher 
tenqierature  at  the  inner  peak  location. 

The  numerical  model  assumes  an  infinitely  fast 
chemical  reaction  rate,  within  the  context  the  flame 
sheet  model.  A  finite  difference  scheme 
(QUICKEST)  is  used  to  obtain  solutions  to  the  2-D. 
time  depoident  Navier-Stokes  equations,  coiqiled  with 
the  continuity  and  mixture  firaction  equations.  The 
simulation  is  seen  to  rqiioduce  the  fcmnation  of  the 
primary  and  secondary  vortex  pair,  and  the  fonnatioa 
of  the  flame  bulge  due  to  the  vortex-flame  interaction. 
The  axial  height  of  the  flame  bulge  is  also  adequately 
reproduced.  Due  to  the  infinitely  fast  chemistry 
assumption,  the  local  quenching  of  the  flame  is  not 
predicted.  The  simulation  nevertheless  shows  a 
double  peak  temperature  profile  at  the  axid  location 
where  local  quenching  was  experimentally  observed. 
At  this  location,  the  gradients  of  temperature  and 
concentration  are  much  sharpo'  than  the  siagle  peak 
profiles.  Due  to  this  sharp  gradient,  the  present 
numerical  grid  resolution  is  not  sufficient  to  precisely 


locate  the  stoichiomeiry.  The  maximum  temperaiure. 
on  contrary  to  the  TFP  data,  is  located  at  the  outer 
peak.  The  predicted  inner  peak  lemperaoire.  however, 
is  much  lower  than  that  observed  in  the  experimenL 

While  the  simulation  shows  the  enhiuiced  mass 
transport  lo  the  reaction  zone  and  the  iianspon 
process  being  important  during  the  vortex  and  flame 
interaction,  it  fails  to  predict  the  flame  stretch  that 
leads  to  local  extinction.  The  much  higher  itmer  peak 
temperature  of  the  TFP  dsu  suggest  that  the  ^le 
rate  chemistry  might  be  important  within  the  vortex, 
a  result  of  the  premixing  of  fuel-product  mixture  with 
the  oxygen  that  leaks  through  the  quenching  location 
or  highly  stretched  the  flame.  Thus,  it  is  possible  to 
observe  distributed  reaction  zones  or  convoluted 
stoichiometric  surfaces  inside  the  ring  vortex  pair. 

This  joint  experimental  and  numerical 
investigation  presenu  an  experimental  system  for 
study  the  vortex-flame  interaction  that  undergoes  a 
thin  flame  regime,  through  a  stretched  flame  regime, 
to  local  exttnctiaii.  The  experimental  results  suggest 
the  possibility  of  partial  premixing  inside  the  ring 
vortex  pair  which  is  yet  to  be  confirmed. 

Several  improvements  of  the  current 
capabilities  are  obvious.  The  finite  rate  chemistry 
must  be  incorporated  in  the  numerical  model  to 
investigare  the  reaction  zone  structure  inside  the 
vortex  ring.  Similarly,  the  experiment  call  for 
species  concentration  roeasutemeats  of  the  ring-vortex 
pair  and  qnaotificatioa  of  die  TFP  time  response  to 
advance  our  undersiandiiig  of  the  type  of  vortex-flame 
interaction  studied  in  this  paper. 
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Figure  1.  Schematics  Driven  Difitekm 
Flame 


Hgure2.  Eiqierimental  Setup  of  Phase- 
Locked  Joint  RMS  and  OH-PUF  Imaging 


Figures.  Experimental  Setup  of 
Thin-Filament  Pyrometiy  Measurement 


Time  (ms) 


Hgure  4.  Hme  History  of  Axial  Velocity  Near 
the  Centerline  of  the  Nozzle  Exit 
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ABSTRACT 

The  crossing  frequency  of  the  inner  vonex  in  a  transitional  jet  diffusion  flame  is  measured 
as  a  function  of  axial  location  using  a  line-visualization  technique.  The  crossing  frequency  is  the 
frequency  for  which  the  vortices  are  convected  across  a  measurement  line  in  a  given  height  in  the 
flame.  Step-wise  decreases  of  the  crossing  frequency  are  observed  in  the  experiments,  in 
qualitative  agreement  with  that  reported  in  the  literature.  The  numerical  simulation  based  on  a 
simple  flame-sheet  model  is  shown  to  reproduce  the  inner  vortex  in  the  transitional  jet  difliision 
flame  studied.  Random  perturbations  applied  at  the  fuel  jet  exit  are  required  for  the  prediction  of 
the  roll-up  vortices  inside  the  flame.  The  simulation  results  also  reveal  the  step-wise  decrease  of 
the  crossing  frequency  along  the  axial  direction.  Hie  step  decrease  of  the  crossing  frequency  is 
shown  to  be  a  result  of  the  vonex  merging  in  jet  diffusion  flames.  The  simulation  captures  the 
major  components  of  the  crossing  frequency  in  the  near  field;  however,  it  is  only  in  qualitative 
agreement  with  the  experimental  crossing  frequency  data  as  a  functitm  of  axial  location.  The 
Strouhal  number  scaling  of  the  vortex  crossing  frequency  shows  that  the  Strouhal  number  is 
bounded  0.2  and  1.0,  typical  of  the  Kelvin-Helmholtz  instability  of  ctrid  jets.  The  present  study 
suggests  that  the  inner  vortex  is  a  Kelvin-Helmholtz  type  of  the  convective  instability  which  is 
modified  due  to  the  presence  of  the  high-temperature  flame  zone.  The  two-dimensional  numerical 
model,  dierefore,  is  capable  of  reproducing  the  frequency  characteristics  of  the  buoyancy-induced 
vortex  in  transitional  jer  diffusion  flames  as  well  as  capturing  the  frequency  components  of  a 
Kelvin-Helmh<dtz  type  of  instalnlity  inside  the  flame. 
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INTRODUCTION 

Two  types  of  flow  vortices  are  known  to  exist  in  transitional  jet  diffusion  flames,  the  large 

toroidal  vortex  outside  the  stoichiometric  surface  (i.e..  the  outer  vonex)  and  the  small  roil-up 

vortex  inside  (i.e..  the  inner  vortex). These  vortices  have  been  illustrated,  for  example,  in  our 

earlier  flow/flame  visualizations^'^  using  a  light  scattering  technique  known  as  the  reactive  Mie 

scattering  (RMS).  The  visualization  shows  that  the  inner  and  outer  vortices  have  different  length 

scales,  for  example,  centimeters  for  the  outer  vortex  and  millimeters  for  the  inner  vwtex.  The 

charactmstic  f  encies  of  these  vortical  structures  are  also  different.  The  outer  vortices  are 

correlated  with  tlame  flicker  frequency typically  in  the  range  10  to  20  Hz.  The  flame  flicker 

is  a  consequence  of  the  buoyancy  acceleration  as  concluded  in  the  numerical  simulations  of  Ellzey 

and  Oran'^  and  Davis  et  al.^  The  microgravity  combustion  experiments  of  Bahadori  et  al.^  also 

reported  a  stationary  (non-oscillating)  flame  in  reduced  gravity.  The  buoyancy  effects  on  flicker 

frequency  have  been  quantified,  for  example,  by  Katta  and  Roquemore^O  and  Sui,^  ^  showing  that 

the  flicker  frequency  is  proportional  to  the  square  root  of  gravitational  acceleradon.  This 

observation  was  in  agreement  with  data  conqiUed  by  Hamins  et  aL^^  in  which  the  Stroohal  number 
vd 

(St  =  d  the  injector  diameter,  v  the  frequency,  u  the  jet  axial  velocity)  was  found  to  be 

prcqxwtional  to  the  square  toot  d  the  Richardson  number  (Ri  =  ^  g  the  gravitatitmal  aoceleration). 

u* 

On  the  inner  v(»tex.  Yule  et  al.^  quantified  the  crossing  frequency  in  the  near-injector 
region  based  on  the  high-speed  Schlieren  movies  taken  for  non-premixed  and  partially-premixed 
flames  of  transitional  jets.  Their  data  showed  that  the  crossing  frequencies  were  typical  about 
several  hundred  hertz  and  the  frequency  value  differed  from  cold  jets  at  a  comparable  Red  (Red  = 
pud/)!;  p  the  density  and  p  the  viscosity,  pn^)eities  based  on  jet-fluid  at  298  K)  condition.  A  step¬ 
wise  decrease  in  vortex-crossing  frequency  along  the  axial  direction  was  also  observed.  The 
coalescence  and  growth  of  inner  vortices  were  identified  to  be  responsible  for  the  step-wise 
decrease  in  crossing  frequency.  In  numerical  simulaticms  of  the  transitional  jet  diffusion  flame,  the 
inner  and  outer  vortices  were  reproduced  by  Ellzey  and  Oran^  and  Katta  and  Roquemore.^^  Katta 
and  Roquemore  also  illustrated  that  the  buoyancy  acceleration  was  important  in  the  maintaining  of 
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the  coherence  of  the  inner  vortices.  While  the  flame  flicker  frequency  or  the  outer  vortex  has  been 
characterized  in  the  literature,  the  experimental  data  to  quantify  the  inner  vortex  in  transitional  jet 
diffusion  flames,  however,  are  lacking.  For  example,  comparison  of  the  simulation  with  the 
experimental  inner-vwtex  crossing  frequency  has  not  been  reported.  The  objectives  of  this  paper 
are  (1)  to  extend  the  data  base  of  the  inner- vortex  crossing  frequency  of  a  non-luminous  diffusion 
flame  and  (2)  to  extend  our  earlier  effort^  of  the  time-accurate,  two-dimensional  simulation  of 
transitional  jet  diffusion  flames  to  compare  the  simulation  with  the  experimental  vortex  crossing 
frequency. 

EXPERIMENTAL  METHODS 

A  nonluminous  jet  diffusion  flame  was  considered  for  the  vortex-crossing  frequency 
measurements.  A  H2/N2  mixture  with  a  molar  ratio  of  3.S  (molecular  weight  of  7.8)  was  used  as 
fuel.  The  stoichiometric  mixture  fraction  of  this  fuel  is  0.126.  The  experiments  were  conducted  in 
an  unconfined,  vertical  combustion  tunnel  in  die  Air  Fence  Wti^t  Laboratory.  This  {qiparatus  has 
been  npotted  elsewhere;^*^  thereftne,  only  a  brief  description  is  given  here.  The  burner  section 
consisted  of  a  vertically  mounted  fuel  tube  surrounded  by  an  aimular-air  jet.  The  annular-air  jet 
provided  a  low-speed  coflow  air,  0.16  m/s  at  the  jet  exit,  to  shield  the  jet  diffusion  flame  from 
ambient  disturbances.  Two  fuel  nozzles  of  5  and  10  mm  in  diameter  were  used.  ITiese  fuel 
nozzles  were  made  from  a  25  mm  stainless-steel  tubing  about  1  m  long,  with  a  t^iered  section 
(length  about  100  mm)  and  a  thin-wall  nozzle  lip  (thickness  less  than  1  mm).  Two  fuel-jet  exit 
velocities  were  examined,  7.8  and  23.0  m/s  for  d  =  10  and  S  mm,  respectively.  The 
corresponding  Reynolds  numbers  were  2300  and  3400. 

The  RMS  visualization  technique  described  in4-6.13  was  employed  for  planar  and  line 
visualizations.  The  planar  visualization  was  for  qualitative  imaging  of  the  vortex  structures  and  the 
line  visualization^^  was  for  quantitative  crossing  frequency  measurements.  A  pulsed  Nd:YAG 
laser  (with  the  532  nm  output)  was  used  as  the  light  source  for  planar  visualization.  The  image 
was  recorded  by  an  intensified  CCD  camera  positioned  at  a  right  angle  from  the  illuminating  laser 
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sheet.  The  line  visualization  technique  was  described  in.  *3  A  continuous-wave  laser  (an  Ar*"  laser 
with  514.5  nm  output)  was  positioned,  passing  through  the  center  of  the  jet  flame  at  axial  heights 
of  z/d  s  0.41  to  6.24  (10  mm  nozzle)  or  2.8  to  12.5  (5  mm  nozzle)  where  z  is  the  axial  distance 
above  the  injector.  A  ten-facet  rotating  mirnM'  was  used  to  focus  to  a  photomultiplier  the  light 
scattering  from  RMS  particles  (Ti02)  along  the  laser  beam,  i.e.,  effectively  scanning  through  the 
laser  beam  over  the  measuring  "volume."  The  photomultiplier  output  was  digitized  by  a  high¬ 
speed  digitizer  (LeCroy  TR8818)  into  256  or  512  data  blocks,  i.e.,  256  or  512  pixels  per  scan. 
The  digitizer  was  operated  at  a  speed  in  the  range  4  to  6  MHz.  The  resultant  spatial  resolution  for 
each  pixel  was  better  than  200  pm.  A  He-Ne  laser  and  a  photodiode  were  used  to  trigger  the 
digitizer  during  each  scan  of  the  minor  facet  Tte  mirror  rotational  q)eed  was  set  in  the  range  200 
to  500  Hz,  corresponding  to  an  effective  sweeping  rate  of  2  to  5  kHz.  This  sweeping  (sampling) 
rate  was  adequate  to  quantify  the  vortex  crossing  frequency  in  this  p^rer.  A  sampling  period  in  the 
range  250  to  500  ms  was  taken  for  FFT  (fast  FtMirier  transform)  analysis  of  the  vortex  crossing 
frequency. 

NUMERICAL  METHODS 

The  explicit  numerical  scheme  of  Wu  and  Chen*^  was  employed.  This  computer  code  is 
similar  to  that  of  Davis  et  al.^,  except  that  a  finite-volume  formulation  was  used  to  replace  the 
QUICKEST  scheme  and  a  different  direct  solver  was  used  for  the  pressure  Poisson  equation.  The 
coiiq>uter  program  solves  the  time-dependent,  axisymmetric  Navier-Stokes  equations  coupled  with 
the  continuity  and  mixture  fraction  equations.  The  mathematical  formulation  is  based  on  a  flame 
sheet  model.  The  major  assumptions  are  an  adiabatic  flame,  a  one-step  irreversible  chemical 
reaction  with  H2O  being  the  only  combustion  product,  an  infrnitely-large  Damkohler  number  or  an 
infinitely-thin  reaction  zone,  equal  diffusion  coefficients  for  all  species,  a  unity  Lewis  number,  and 
a  low  Mach  number,  near-laminar,  Newtonian  flow  with  negligible  radiative  heat  transfer  and 
viscous  dissipation.  The  present  model  is  similar  to  that  employed  by  Katta  et  al.^5  in  the 
simulation  of  a  similar  H2/N2  diffusion  flame  studied  in  this  paper.  The  governing  equations  and 
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the  boundary  conditions  will  not  be  repeated  here  because  they  have  been  documented^*^,  and  in 
detail  by  Lee.'^ 

Hie  near-field  development  of  a  jet  flow  is  known  to  be  sensitive  to  the  jet  initial  shear 
layer.  An  initial  momentum  thickness  of  0.21  mm  using  a  hyperbolic-tangent  profile  was 
specified.  This  momentum  thickness  was  similar  to  that  employed  by  Grinstein  et  al.,^^  and 
consistent  with  the  measured  thickness  0.3  mm  determined  from  the  LDA  data.  The  same  shape 
of  the  hyperbolic-tangent  profile  was  ^lecified  for  the  mixture  fraction  at  the  jet  exit,  with  pure  fuel 
having  a  value  1  and  pure  air  a  value  0.  Small  perturbation  imposed  to  the  axial  velocity  at  the  jet 
exit  was  necessary  to  initiate  the  roll-up  of  the  inner  vortices  in  the  computation.  A  random 
number  generator  was  used  to  introduce  white  noise  at  the  fuel-jet  exit,  over  a  region  of  4  to  8  grid 
points  within  2  mm  in  the  jet  shear  layer.  The  perturbation  an^litude  was  set  not  to  exceed  4  %  of 
local  velocity.  The  amplitude  and  the  grids  were  randomly  selected. 

To  resolve  the  small  length  scale  associated  with  the  irmer  vortex,  180  grid  points  were 
used  for  an  axial  distance  of  600  mm  and  80  grid  points  for  a  radial  distance  of  100  mm.  The  grids 
in  the  radial  direction  woe  concentrated  in  the  shear  layer  region  and  around  the  stoichiometric 
surface.  The  grid  spacing  in  axial  directirm  was  gradually  increased  as  z  was  increased.  The 
boundaries  were  set  far  from  the  stoichkxnetiic  surface  to  prevent  the  propagation  of  disturbance 
due  to  inaccuracies  in  specifying  the  boundary  conditions  of  the  finite  computational  domain.  The 
time  step  was  set  to  satisfy  the  Courant  condition.  Typically,  a  time  step  of  15  ps  or  shorter  was 
used.  This  time  step  was  sufficient  to  reserve  the  evolution  of  the  inner  vortex.  The  computation 
was  executed  for  25,000  time  steps  for  each  condition  considered.  A  computer  workstation 
(Apollo  DN  KXXX))  was  used,  15  CPU  hours  was  required  for  25,(X)0  time  steps  (or  375  ms  of 
real  time). 

RESULTS  AND  DISCUSSION 

Hie  present  injector  yields  "top-hat"  velocity  profiles  at  the  burner  exit  as  quantified  by  the 
LDA  measurements.^®  These  "top-hat"  velocity  profiles  can  be  approximated  by  the  hyperbolic- 
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tangent  profiles  used  in  the  simulation.  A  typical  RMS  visualization  of  the  flames  studied  is 
illustrated  in  Hg.  1.  Both  the  outer  and  inner  vortices  can  be  seen  in  this  visualization.  The  laser 
beam  for  line  visualization  was  also  illustrated  (not  scaled)  in  this  figure.  The  line  visualization  of 
a  similar  flame  with  d  =  10  mm  is  shown  in  Fig.  2.  Two  flow  structures  were  identified:  outer  and 
inner  vortices,  marked  by  scattered  light  from  Ti02  particles  outside  or  inside  the  stoichiometric 
surface.  The  outer  vortex  exhibits  a  crossing  frequency  of  15.6  Hz  from  the  FFT  analysis,  cf. 
Fig.  3,  in  agreement  with  the  IS  Hz  flame-flicker  frequency  obtained  for  this  flame  using  a  laser- 
deflection  technique,  The  small  cells  are  due  to  the  crossings  of  the  inner  vonex  visualized  in 

the  space-time  coordinate  of  Fig.  2.  These  cells  registered  the  signals  from  the  discrete  and  finite¬ 
sized  vortex  when  passing  throu^  a  fixed  line  in  space.  A  frequency  around  420  Hz  was  obtained 
from  the  PSD  (power  spectral  density)  which  was  shown  in  Fig.  3.  This  crossing  frequency  was 
independent  of  the  spatial  location  within  the  vortex  cells.  In  general,  interrogation  over  25  pixels 
(or  about  2  to  3  mm  along  the  laser  beam)  was  made  to  obtain  the  spectrum  in  the  present  study. 

To  assess  whether  the  numerical  simulation  can  reproduce  the  crossing  frequency  of  the 
inner  vortex,  the  line  visualization  of  the  sunulation  is  shown  in  Hg.  4  for  the  condition  of  Fig.  2, 
using  the  predicted  radial  velocity.  The  fine  visualization  of  the  Emulation  showed  the  crossings  of 
the  inner  vortex,  marked  by  the  cells  in  die  time-space  comditutie.  This  visualizatimi  was  similar  to 
the  experimental  observation.  The  simulation  also  showed  that  die  vortex  crossing  frequency  was 
independent  of  ^atial  locations  when  the  interrogation  was  taken  within  the  cells.  Specifically  for 
d  =  10  mm,  the  predominant  frequency  was  nearly  identical  for  0.29  <  r/d  <  0.71  (r  the  radial 
distance),  consistent  with  the  experimental  observation.  A  standard  location  at  r  =  3.98  mm  was 
chosen  for  analysis  of  the  simulation  results.  The  line  visualization  also  showed  that  the  vortex  cell 
grows  in  size  at  downstream  locations.  At  downstream  locations,  the  crossing  frequency  seems  to 
have  decreased. 

To  compare  the  PSD  at  different  axial  locations.  Fig.  5  summarized  the  experimental  and 
simulation  results  at  three  axial  locations.  Different  vortex  crossing  frequencies  were  found  in  the 
spectrum.  For  example,  the  measurements  re  pealed  two  major  components  centered  around  460 
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Hz  and  840  Hz  at  z/d  =  2.0,  one  component  around  490  Hz  at  z/d  =  4.0  and  320  Hz  at  z/d  =  5. 1. 
At  z/d  =  4.0  and  5.1,  the  flicker  frequency  around  15  Hz  was  also  registered.  The  simulation 
showed  that  the  vortex  crossing  frequency  peaked  around  400  Hz  at  z/d  =  2.0  and  300  Hz  at  z/d  = 
3.6  and  z/d  =  4.9.  The  flame  flicker  frequency  around  15  Hz  was  also  predicted  in  the  simulation. 
These  PSD's  showed  that  the  vortex  crossing  frequency  decreases  as  axial  distance  is  increased. 

To  assess  whether  the  step  change  described  by  Yule  et  al.^  also  exists  in  the  flame  studied, 
the  inner  vortex  crossing  frequency  was  plotted  as  a  function  of  z/d  in  Fig.  6.  The  strongest 
component  in  the  spectrum  other  than  the  flicker  frequency  was  taken  as  the  crossing  frequency. 
Indeed,  step  change  was  evident  For  example,  the  measured  frequency  decreased  from  1400  Hz 
to  850  Hz  at  z/d  =  1,  to  450  Hz  at  z/d  =  2,  and  to  300  Hz  at  z/d  -  5.  At  z/d  =  6,  there  seems  to  be 
another  step  change  of  the  frequency  from  300  Hz  to  250  Hz,  although  this  change  is  not  as 
evident  as  the  above-noted  decreases.  The  numerical  simulation  also  predicted  a  similar  step 
decrease  in  the  crossing  frequency.  The  simulation  conastently  underestimated  the  vertex  crossing 
frequency,  although  the  discrepancy  was  reduced  at  downstream  locations.  The  comparison  was 
not  extended  to  downstream  locations  because  the  inner  vtxtex  has  beccxne  three-dimensional  at  z/d 
around  8.  The  results  of  the  vortex  crossing  frequency  in  the  5  mm  diffusion  flame  are  similar  to 
the  10  mm  flame  and  step-wise  decreases  are  also  obsoved.^^ 

To  assist  the  interpretation  of  step  decrease  of  the  crossing  frequency,  the  predicted  time 
trace  of  radial  velocity  at  four  axial  heights  is  shown  in  Hg.  7  for  the  10  mm  flame.  The  velocity  at 
the  near-injector  exit  (z/d  =  0.93)  appeared  to  be  noisy;  the  periodicity  in  the  signal  was  not 
obvious.  A  "modulation"  or  "merging"  of  the  signal  into  a  more  periodic  time  trace  seems  to  have 
occurred  at  downstream  locations.  For  example,  the  time  trace  merged  into  six  cycles  at  z/d  =  0.93 
over  the  time  period  9  to  12  ms  (cf..  Fig.  7),  into  two  cycles  at  z/d  =  1.54  over  10  to  13  ms,  and 
into  a  single  cycle  at  z/d  =  3. 12  over  13  to  16  ms.  The  time-trace  results  suggested  that  the  vortex 
merging  was  responsible  for  the  step  decrease  of  the  crossing  frequency.  The  vortex  merging  in 
transitional  jet  diffusion  flames  has  been  identifred,  for  example,  in  RMS  visualizations.  ^'4  The 
vortex  merging  or  step-wise  frequency  decrease  in  jet  diffusion  flames  is  similar  to  the  vortex 
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pairing  (or  period  doubling)  of  the  cold  jets.^l  The  merging  process  in  jet  diffusion  flames, 
however,  is  different  from  cold  jets.  A  step-wise  decrease,  not  equal  to  the  frequency  halving  of 
the  cold  jets,  was  typical  in  jet  diffusion  flames.  The  jet  diffusion  flame  also  has  a  complex 
spectrum  showing  a  txoad-band  around  the  predominant  frequency  (cf..  Fig.  5)  while  a  simple  and 
well  defined  peak  is  typical  for  cold  jets.  The  broadening  of  the  spectrum  is  due  to  the  presence 
of  the  high-temperature  flame  zone  which  limits  the  growth  of  the  inner  vortex.  The  vortex¬ 
crossing  frequency  is  further  complicated  by  the  time-varying  flame  zone  due  to  the  buoyant 
instability  (and  other  sources)  and  the  accompanied  time-varying  velocity.  As  a  result,  a  more 
complex  spectrum  is  typical  for  the  jet  diffusion  flame  compared  to  the  cold  jet  at  a  similar  Re^ 
cofKlition. 

The  Strouhal  number  has  been  used  in  the  scaling  of  the  Kelvin-Helmholtz  instability  in 
non-reacting  jets.  The  Strouhal  based  on  the  injector  diameter  (Std)  was  plotted  as  a  function 

of  Rcz  (based  on  axial  distance)  in  Fig.  8,  including  three  more  data  sets  of  simulation  results.  All 
the  conditions  in  Fig.  8  have  a  Red  in  the  range  2300  to  3400.  The  results  showed  that  the  Std 
generally  in  the  range  0.2  to  1.0,  and  Std  decreased  as  Rez  was  increased.  A  nearly  constant  value 
around  0.2  was  found  for  Rcz  in  the  range  50,000  to  60,000.  The  Std  of  0.2  was  typical  for  a 
non-reactive  jet  with  a  thin  initial  shear  layer  (or  a  top-hat  velocity  profile).  When  the  momennim 
thickness  (6)  was  used  as  the  length  scale,  Ste  was  shown  to  be  bounded  by  0.005  and  0.05  for 
the  conditions  examined.  These  observations,  along  with  the  fact  that  a  random  perturbation  was 
needed  for  simulation  of  the  inner  vortex,  the  inner  vortex  appears  to  be  a  Kelvin-Helmholtz  type 
of  the  convective  instability  modified  due  to  the  presence  of  a  high  temperature  flame  zone. 

SUMMARY  AND  CONCLUSIONS 

The  crossing  frequency  of  the  inner  vortex  in  a  transitional  jet  diffusion  flame  (with  a  near 
"top-hat"  velocity  prc^e  at  the  burner  exit)  is  studied.  The  crossing  frequency  has  been  measured 
for  two  nitrogen-diluted  hydrogen  jet  diffusion  flames,  employing  a  line  visualization  technique. 
Spatially-resolved  vortex  crossing  frequencies  are  documented.  Step-wise  decreases  of  the 
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crossing  frequency  are  observed  for  the  conditions  examined,  in  qualitative  agreement  with  Yule  et 
al.^  The  numerical  simulation  based  on  a  simple  flame-sheet  model  is  shown  to  be  capable  of 
reproducing  the  outer  and  the  inner  vortices  in  the  flame  studied.  Random  perturbation  imposed  to 
the  axial  velocity,  however,  is  needed  at  the  fuel  jet  exit  to  initial  the  roll-up  of  the  inner  vortices. 
The  simulation  also  reveals  step-wise  decreases  of  the  crossing  frequency  along  the  axial  direction. 
The  step  decrease  of  the  crossing  frequency  has  been  shown  to  be  a  result  of  the  vortex  merging  in 
jet  diffusion  flames.  The  broad-band  spectrum  is  believed  to  be  due  to  the  presence  of  the  flame 
zone.  The  predicted  crossing  frequencies,  although  are  only  in  qualitative  agreement  with  the 
spadally-specific  experimental  data,  do  capture  the  majOT  components  in  the  near  field  of  the  flame. 

The  Strouhal  number  scaling  of  the  vortex  crossing  frequency  (based  on  the  jet  diameter) 
versus  the  Reynolds  number  based  on  axial  distance  shows  that  the  Strouhal  number  is  bounded 
by  0.2  and  1.0,  typical  of  the  Kelvin-Helmholtz  instability  of  cold  jets.  This  fact  along  with 
simulation  requiring  a  perturbation  source  at  the  jet  exit  (which  was  not  needed  for  the  buoyant 
instability^)  suggests  that  the  inner  vortex  is  a  Kelvin-Helmholtz  type  of  the  ccmvective  instability, 
and  the  instability  was  modified  due  to  the  presence  of  the  high-temporature  flame  zone.  Further 
validation  of  this  hypothesis,  however,  is  needed.  Hie  present  study,  in  our  opinion,  represents  a 
significant  advancement  in  our  effort  to  verify  the  two-dimensional,  time-accurate  prediction  of 
transitional  jet  diffusion  flames.  The  numerical  model  is  shown  not  only  it  can  reproduce  the 
frequency  characteristics  of  the  buoyancy-induced  vortex  in  transitional  jet  diffusion  flames®  but 
also  it  can  capture  the  frequency  components  of  a  Kelvin-Helmholtz  type  of  instability  inside  the 
flame. 

There  are  several  limitations  in  the  present  numerical  model,  for  example,  the  unity  Lewis 
number  assumption  and  the  two-dimensional  simulation.  The  non-unity  Lewis  number  effect  has 
been  shown  to  alter  the  local  temperature  due  to  the  preferential  diffusion  in  the  flame  studied,^^ 
but  it  did  not  alter  the  global  dynamics  of  the  flame.  To  what  extent  the  preferential  diffusion 
affects  the  time  characteristics  of  the  inner  vortex  is  yet  to  be  ascertained.  To  apply  the  time- 
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accurate  calculation  for  highly  three-dimensional  turbulent  flames,  of  course,  remains  a  continued 
challenge  to  the  authors. 
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LIST  OF  FIGURES 


Fig.  1.  RMS  visualization  of  a  H2/N2  jet  diffusion  flame  (uj  =  23.5  m/s,  d  =  5  mm,  Rcj  =  3400; 
Ua  =  0.16  m/s);  Uj  the  fuel-jet  exit  velocity,  d  the  diameter,  Rej  the  jet  exit  Reynolds 
number  and  Ua  the  coflow  annular  air  velocity  at  the  jet  exit;  the  image  view  was  about 
65  mm  (radial)  by  200  mm  (axial);  the  high  temperature  surface  was  inferred  from  thin- 
filament  visualization  measurements.^^ 

Fig.  2.  Line  visualization  of  a  H2/N2  jet  diffusion  flame  (uj  =  7.8  m/s,  d  =  10  nun,  Rej  =  2300; 
Ua  =  0.16  m/s);  the  laser  beam  was  located  at  z/d  =  4.7,  and  the  minor  scanning  rate  was 
set  at  4  kHz  and  the  digitizing  rate  at  6  MHz. 

Fig.  3.  The  PSD  of  data  shown  in  Hg.  2. 

Fig.  4.  line  visualization  of  the  numerical  simulation  a  H2/N2  jet  diffusion  flame  (uj  =  7.8  m/s, 
d  =  10  mm.  Red  =  23(X);  Ua  =  0.16  m/s). 

Fig.  5.  The  PSD  ai  a  H2/N2  jet  diffusion  flame  (uj  =  7.8  m/s,  d  =  10  mm.  Red  =  2300;  Ua  =  0.16 
m/s)  due  to  the  inner  vortex  crossing  (a)  experimental  data  at  zjd  -  2.0, 4.0  and  5.1.,  and 
(b)  simulation  results  at  z/d  =  2.0, 3.6  and  4.9. 

Fig.  6.  The  vortex  crossing  frequency  of  a  H2/N2  jet  diffusion  flame  (uj  =  7.8  m/s,  d  =  10  mm. 
Red  =  2300;  Ua  =  0.16  m/s)  as  a  function  of  z/d;  experimental  data  shown  by  closed 
symbols,  simulation  by  open  symbols. 

Fig.  7.  The  time  trace  of  radial  velocity  at  different  axial  locationsf  a  H2/N2  jet  diffusion  flame 
(uj  =  7.8  m/s,  d  =  10  mm.  Red  =  2300;  Ua  =  0.16  m/s);  the  circled  regions  illustrate  the 
evolving  (or  noodulation)  of  the  time  trace  into  periodical  signals. 

Fig.  8.  Scaling  of  the  crossing  frequency  by  the  Strouhal  number  based  on  nozzle  diameter,  Std, 
versus  the  Reynolds  number  based  on  axial  location  Rcz,  for  different  flow  conditions 
and  different  axial  locations;  *  denotes  experimental  measurements. 
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